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ABSTRACT 

We present tiie galaxy stellar mass function (GSMF) to redshift z ~ 1, based on the analysis of about 8500 galaxies with / < 22.5 
(AB mag) over 1.4 deg^, which are part of the zCOSMOS-bright 10k spectroscopic sample. We investigate the total GSMF, as 
well as the contributions of early- and late-type galaxies (ETGs and LTGs, respectively), defined by different criteria (broad- 
band spectral energy distribution, morphology, spectral properties, or star formation activities). We unveil a galaxy bimodality 
in the global GSMF, whose shape is more accurately represented by 2 Schechter functions, one linked to the ETG and the other 
to the LTG populations. For the global population, we confirm a mass-dependent evolution ("mass-assembly downsizing"), i.e., 
galaxy number density increases with cosmic time by a factor of two between z — 1 and z = for intermediate-to-low mass 
{\og{M/MQ)^ 10.5) galaxies but less than 15% for \og{M/MQ)> 11. We find that the GSMF evolution at intermediate-to- 
low values of A4 {\og{A4 / A4q)< 10.6) is mostly explained by the growth in stellar mass driven by smoothly decreasing star 
formation activities, despite the redder colours predicted in particular at low redshift. The low residual evolution is consistent, 
on average, with ~ 0.16 merger per galaxy per Gyr (of which fewer than 0.1 are major), with a hint of a decrease with cosmic 
time but not a clear dependence on the mass. From the analysis of different galaxy types, we find that ETGs, regardless of the 
classification method, increase in number density with cosmic time more rapidly with decreasing A4, i.e., follow a top-down building 
history, with a median "building redshift" increasing with mass (z > 1 for log(A1/A4o)> 11), in contrast to hierarchical model 
predictions. For LTGs, we find that the number density of blue or spiral galaxies with \og{M/AiQ)> 10 remains almost constant 
with cosmic time from z ^ 1. Instead, the most extreme population of star-forming galaxies (with high specific star formation), 
at intermediate/high-mass, rapidly decreases in number density with cosmic time. Our data can be interpreted as a combination 
of different effects. Firstly, we suggest a transformation, driven mainly by SFH, from blue, active, spiral galaxies of intermediate 
mass to blue quiescent and subsequently (1-2 Gyr after) red, passive types of low specific star formation. We find an indication 
that the complete morphological transformation, probably driven by dynamical processes, into red spheroidal galaxies, occurred on 
longer timescales or followed after 1-2 Gyr. A contiimous replacement of blue galaxies is expected to be accomplished by low-mass 
active spirals increasing their stellar mass. We estimate the growth rate in number and mass density of the red galaxies at different 
redshifts and masses. The corresponding fraction of blue galaxies that, at any given time, is transforming into red galaxies per 
Gyr, due to the quenching of their SFR, is on average ~ 25% for \og{M/MQ)< 11. We conclude that the build-up of galaxies and 
in particular of ETGs follows the same downsizing trend with mass (i.e. occurs earlier for high-mass galaxies) as the formation of 
their stars and follows the converse of the trend predicted by current SAMs. In this scenario, we expect there to be a negligible 
evolution of the galaxy baryonic mass function ( GBMF) for the global population at all masses and a decrease with cosmic time 
in the GBMF for the blue galaxy population at intermediate-high masses. 
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Tracing the history of both galaxy star formation (Lilly 
et al. 1996; Madau et al. 1996; Madau, Pozzetti, Dickinson 
1998) and stellar mass assembly (Dickinson et al. 2003) over 
cosmic time represent major challenges in modern cosmol- 
ogy. In particular, it is still uncertain how the bimodality 
seen in the local Universe in terms of galaxy properties (e.g., 
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colour, morphology, star formation, and spectral features; 
KaufTmann et al. 2003, Baldry et al. 2004, Brinchmann et 
al. 2004) evolves and when it was created. Galaxies in the 
local Universe exhibit distinctive bimodal colour distribu- 
tions (Strateva et al. 2001, Hogg et al. 2002, Blanton et al. 
2003), which suggests that different evolution histories exist 
for galaxies lying on the two sequences (Menci et al. 2005; 
Scarlata et al. 2007b; De Lucia et al. 2007). Even if colour 
bimodality has already been observed and studied at higher 
redshift (BeU et al. 2004, Weiner et al. 2005, up to z - 1; 
Franzetti et al. 2007, Cirasuolo et al. 2007, up to z ~ 1.5; 
Giallongo et al. 2005, Cassata et al. 2008, Williams et al. 
2009, up to 2; ~ 2), no study until now has fully evaluated 
how the number densities of the two populations compare 
and evolve with cosmic time. 

It is still unclear, indeed, how the galaxy stellar mass 
function (GSMF) is linked to the galaxy bimodality. Baldry 
et al. (2004, 2006, 2008) noted a bimodal shape in the lo- 
cal GSMF in the SDSS with an upturn at masses lower 
than 10^ Mq, related to the two different galaxy popula- 
tions. It is, therefore, extremely interesting to explore how 
this shape evolves with redshift and which mechanisms con- 
tribute to its appearance. From previous surveys, it is well 
established that there has been an increase/decrease in the 
fraction of red/blue or early/late-type galaxies with cosmic 
time at intermediate masses since z ~ 1 (Fontana et al. 
2004; Bundy et al. 2006; Arnouts et al. 2007; Scarlata et al. 
2007b; Vergani et al. 2008), while the precise evolution in 
absolute number densities of the two populations remains 
controversial, in particular for intermediate/massive galax- 
ies. A knowledge of this evolution may help to constrain 
the quenching mechanisms responsible for downsizing. 

In agreement with the first formulation introduced by 
Cowie et al. (1996), a downsizing scenario in both age and 
star formation history has been proposed by several ob- 
servational studies, i.e., more massive galaxies form their 
stars earlier and more rapidly than lower mass ones (here- 
after age downsizing: Brinchmann & Ellis 2000; Gavazzi & 
Scodeggio 1996; Fontana et al. 2004, Kodama et al. 2004; 
Bauer et al. 2005; Feulner et al. 2005a, b; Juneau et al. 
2005; Borch et al. 2006; Cucciati et al. 2006; Vergani et 
al. 2008). Age downsizing occurs separately for each of the; 
two populations causing galaxy bimodality, i.e., a red peak 
(Thomas et al. 2005, Fontana et al. 2004, Thomas et al. 
2009) and a blue one (Noeske et al. 2007a, b). However, it 
remains unclear whether the age downsizing is coupled with 
a mass-assembly downsizing scenario for galaxy evolution 
and formation (Fontana et al. 2004, 2006; Pozzetti et al. 
2003, 2007; Cimatti et al. 2006; Bundy et al. 2006), i.e., 
if the more massive galaxies assembled their mass earlier 
than lower mass ones. Furthermore, do low-mass galaxies 
contain younger stars and assemble later even within the 
same spectral type? The hierarchical model of De Lucia et 
al. (2006) predicts, for example, a bottom-up assembly his- 
tory for elliptical galaxies (also called "upsizing") following 
the hierarchical growth of dark matter haloes, in contrast to 
a top-down, downsizing scenario for the formation of their 
stars. 

A fundamental role towards answering these questions 
is played by deep surveys, which sample thousands of galax- 
ies across large portions of the sky. Very deep surveys have 
been exploited to describe the shape of the stellar mass 
function at high redshift (Fontana et al. 2006; Drory et al. 
2005; Gwyn & Hartwick 2005; Bundy et al. 2006; Pozzetti 



et al. 2007), but a clear picture about stellar mass assem- 
bly and how it depends on mass (mass-assembly downsiz- 
ing) and galaxy type has not yet emerged. Previous stud- 
ies have explored the evolution of different galaxy types in 
deep near-IR surveys, such as K20, by means of the _ftr-band 
luminosity function (Pozzetti et al. 2003) and the galaxy 
stellar mass function (Fontana et al. 2004), using spectral 
classification (i.e., absorption-line galaxies versus emission- 
line galaxies), and in larger optical and near-IR surveys 
such as VVDS, C0MB017, and DEEP2, using colours or 
spectra to define galaxy types (Bell et al. 2004; Cimatti et 
al. 2006; Faber et al. 2007; Zucca et al. 2006; Arnouts et 
al. 2007; Vergani et al. 2008). Even if the results of these 
surveys remain disputed (compare for example Bell et al. 
2004 and Cimatti et al. 2006 for the same dataset), most 
of these studies agree that luminous and rather massive 
old galaxies were already quite common at z 1 and that 
their number density declines rapidly at yet higher redshift. 
This suggests that merger events are ruled out as the ma- 
jor mechanism behind their assembly history below z ~ 1. 
However, observational results about major merging and 
dry merging are still contradictory (see Bell et al. 2006; 
van Dokkum 2005; Lin et al. 2004; de Ravel et al. 2009; 
and Renzini 2007 for a summary). It is also known that 
the less luminous/massive ETGs decline in number den- 
sity steadily with redshift (Cimatti et al. 2006). Scarlata et 
al. (2007b) show similar results for the photometric survey 
COSMOS, using both morphologically and photometrically 
selected subsamples of early- type galaxies. The evolution in 
number density of the most massive late-type star-forming 
galaxies continues, however, to be unclear. 

Another open question is whether the evolution of the 
observed GSMF with cosmic time is driven mainly by ei- 
ther merging events at any given mass, as predicted by 
hierarchical galaxy formation models (Cole et al. 2000; 
Menci et al. 2005; Bower et al. 2006; De Lucia et al. 2006; 
Monaco et al. 2006) or star formation histories (SFHs, see 
Vergani et al. 2008; Walcher et al. 2008). It is widely be- 
lieved that galaxies are assembled by hierarchical mergers 
between massive cold dark matter haloes, in which bary- 
onic star-forming matter is embedded. However, most of 
the hierarchical galaxy assembly models are unable to com- 
pletely account for the observed GSMF and its evolution 
(see Fontana et al. 2004, 2006, Monaco et al. 2006, Caputi 
et al. 2006, Marchesini et al. 2009, Fontanot et al. 2009 for a 
detailed comparison with models). For instance, some mod- 
els tend to underpredict the high-mass tail and overpredict 
the rate of its evolution (Fontana et al. 2004, 2006), even 
in the extreme case of evolution driven purely by mergers 
(Monaco et al. 2006). Only the assumption that a signif- 
icant fraction, ~ 30 %, of stars are scattered within the 
diffuse stellar component at each merger event leads to 
significant suppression of the predicted evolution rate, in 
closer agreement with observational constraints (Monaco et 
al. 2006). On the other hand, most of the models overpre- 
dict the number density of relatively low-mass galaxies (see 
Fontana et al. 2006, Kitzbichler & White 2007, Marchesini 
et al. 2009, Fontanot et al. 2009). Furthermore, the models 
do not reproduce the downsizing trend in stellar mass ob- 
served for elliptical galazies (Cimatti et al. 2006). According 
to the latest observational results, the current galaxy for- 
mation and evolutionary scenario is becoming one in which 
a smoother evolution in mass growth and star formation 
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(due to cold gas accretion) plays a major, if not dominant, 
role compared to dark matter (major) merging events. 

In this paper, we use the zCOSMOS spectroscopic sur- 
vey (Lilly et al. 2007, 2009) to complete a comprehensive 
analysis of the GSMF. Compared to previous spectroscopic 
surveys, the larger area and number of spectroscopic red- 
shifts of zCOSMOS allows higher precision, higher quality 
statistics, and lower cosmic variance (Lilly et al. 2009) to 
be achieved in our analysis, in particular for the massive 
end of the GSMFs. We derive the GSMF and its evolution 
with cosmic time since z = 1 , as well as the contribution of 
different galaxy populations, using numerous classification 
methods defined in terms of their colours, morphologies, 
star formation activities, or spectroscopic classifications. 
By using spectroscopic redshifts, we are able to study the 
shape of the global GSMF to high precision, its evolution, 
and how it is related to the bimodalities in galaxy proper- 
ties (colours, morphologies, spectral properties) observed to 
high redshifts. We explore the different roles of merging and 
SFHs as a function of and cosmic time to explain the 
observed GSMF evolution. We also explore the evolution 
with cosmic time of the GSMF for different galaxy types 
and propose an evolutionary scenario. 

Throughout the paper, we adopt the cosmology ilm = 
0.25 and = 0.75, with hjo = Hq/{70 km s'^ Mpc^^). 
Magnitudes are given in the AB system. 



2. COSMOS and the bright 10k zCOSMOS 
spectroscopic sample 

The Cosmic Evolution Survey (COSMOS, Scoville et al. 
2007) is the largest HST survey (640 orbits) ever under- 
taken, imaging a field of '--^ 2 deg^ with the Advanced 
Camera for Surveys (ACS) with F814W (Koekemoer et al. 
2007). This survey was designed to probe galaxy evolution 
and the effects of environment to high redshift. COSMOS 
observations include good coverage of the field with multi- 
band photometry from the UV (with GALEX, Zamojski et 
al. 2007), optical (with Subaru and CFHT, Taniguchi et al. 
2007, Capak et al. 2007), NIR (with CTIO, KPNO, Capak 
et al. 2007, and CFHT, McCracken et al. 2010), to MIR 
and FIR (S-COSMOS with Spitzer, Sanders et al. 2007), 
in combination with a multiwavelength dataset from ra- 
dio (with VLA, Schinnerer et al. 2007), millimetre (with 
MAMBO-2 at the IRAM telescope, Bertoldi et al. 2007), 
to X-rays (with XMM, Hasinger et al. 2007, and Chandra, 
Elvis et al. 2009). 

The zCOSMOS spectroscopic survey (Lilly et al. 2007) 
is an ongoing ESO Large Programme (~ 600 hours of ob- 
servations) aiming to map the COSMOS field with the 
Vi sible M uhi-Object Spectrograph (VIMOS, Le Fevre et 
al. I2003a|) . mounted on the ESO Very Large Telescope 
(VLT). The zCOSMOS survey consists of a bright part, 
with spectroscopy limited to objects in the magnitude range 
15.0 < / < 22.5, and of a deep part, which measures red- 
shifts of i? < 25.25 of galaxies colour-selected to be in the 
range 1.4 < 2; < 3, within the central 1 deg^. The bright 
part has already produced redshifts and spectra for about 
10 000 galaxies over ~ 1.4 deg^, the so-called lOk-bright 
spectroscopic sample (Lilly et al. 2009), with an average 
sampling rate of about 33%. For more details about the 
zCOSMOS lOk-bright sample, we refer to Lilly et al. (2009). 
Here we recall that the VIMOS spectroscopic observations 



were completed using the red medium resolution R ^ 600 
MR grism (5550-9650 A), the spectra were reduced using 
the VIMOS Interactive Pipeline Graphical Interface soft- 
ware (VIPGI, Scodeggio et al. 2005) and redshift measure- 
ments were visually determined after a first estimate had 
been provided by an automatic package (EZ, Garilli et al. 
in preparation) . A quality flag was assigned to each redshift 
measurement. This fiag ranges from (failed measurement) 
to 4 (100% confidence level), and fiag 9 indicates spectra 
with a single emission line, for which multiple redshift so- 
lutions are possible. Additional details about the spectro- 
scopic quality fiags and their probability of being robust are 
given in Lilly et al. (2009; see their Table 1). In addition to 
the confidence classes described in Lilly et al. (2007), a deci- 
mal place (from 5 to 1, see Table 3 in Lilly et al. 2009) in the 
class is added to indicate the level of consistency between 
the spectroscopic and photometric redshifts obtained by the 
Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA; 
Feldmann et al. 2006), using the optical to infrared SED. 

2.1. The selected galaxy spectroscopic sample 

The analysis presented in this paper is based on the 
zCOSMOS-lOk bright sample (Lilly et al. 2009). From the 
total sample of 10 644 objects observed spectroscopically, 
we used the objects within the statistical sample defined 
in the magnitude range 15 < / < 22.5, and removed the 
spectroscopically confirmed stars, broad-line AGNs, and 
the galaxies with low quality redshift fiag (fiag< 1.5, i.e., 
with a verification rate < 90% and a spectroscopic redshift 
inconsistent with the photometric redshift). To achieve re- 
liable SED fitting to multi-band photometric data, we ex- 
cluded objects with apparent magnitudes measured in fewer 
than 5 bands (~ 1.7%) and objects for which the ground 
photometry can be affected by the blending of sources, as 
inferred from the number of ACS sources brighter than 
/ = 22.5 within 0.6" (- 0.5%). This resuhed in the se- 
lection of 8450 galaxy spectra with secure spectroscopic 
measurements (7936 in the redshift range where the fol- 
lowing analysis is carried out, z = 0.1 — 1) over 5045.80 
arcmin^. Objects with redshift fiags < 1.5 are taken into 
account statistically (see Sect.jSJ and Bolzonella et al. 2009 
and Zucca et al. 2009 for details). 

The final selected sample of 8450 galaxies, which have 
highly reliable spectroscopic redshifts, have complete multi- 
band photometric coverage (from UV to IRAC) and mor- 
phological classifications, along with morphological param- 
eters from both the Zurich Estimator of Structural Types 
(ZEST; Scarlata et al. 2007a) and estimates obtained in 
Marseille (Cassata et al. 2008, Tasca et al. 2009, MRS here- 
after) . 

2.2. Photometric data 

The COSMOS field has been covered by multiband pho- 
tometry over a wide range of wavelengths. In this pa- 
per, we used the observed magnitudes in 10 photomet- 
ric bands (CFHT u* and Ks, Subaru Bj, Vj, g+ , r+, 
and z+, and Spitzer IRAC at 3.6 /im and 4.5 fim). 
Descriptions of the photometric catalogs are given in Capak 
et al. (2007), Sanders et al. (2007), and McCracken et al. 
(2010). Following the same approach as Capak et al. (2007; 
see their Table 13), the photometry was adjusted by ap- 
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plying fixed zeropoint offsets to the observed magnitudes 
in eacli band to statistically reduce tlie differences between 
observed and reference magnitudes computed from a set 
of template SEDs, finding in general for each band very 
similar offsets. 

3. Estimate of the stellar masses 

We used the stellar masses {A4) estimated from a fit to the 
multicolour spectral energy distribution (SED), using the 
observed magnitudes in 10 photometric bands from m* to 
4.5 /xm (see Sect. 12. 2) ). and following the method described 
in Pozzetti et al. (2007). In a parallel paper, Bolzonella et al. 
(2009) describe the different methods used to compute stel- 
lar masses, based on different assumptions about the pop- 
ulation synthesis models (Bruzual & Chariot 2003, BC03; 
Maraston 2005, M05; Chariot & Bruzual 2007, CB07) and 
the star- formation histories (SFHs), such as smooth expo- 
nentially decreasing or complex SFHs with the addition of 
secondary bursts. In Bolzonella et al., the associated un- 
certainties and degeneracies are also discussed (see also 
Fontana et al. 2004, Pozzetti et al. 2007, Marchesini et al. 
2009). Here we just recall that the accuracy of the pho- 
tometric stellar masses is satisfactory overall, with typical 
dispersions caused by statistical uncertainties and degen- 
eracies of the order of 0.2 dex. We note that the addition 
of secondary bursts to a continuous star-formation history 
produces systematically higher (up to 40% on average) stel- 
lar masses (Pozzetti et al. 2007), and in particular Fontana 
et al. (2004) demostrated that the effect is larger when the 
fraction of mass produced in the burst is small. On the other 
hand, population synthesis models with a TP-AGB phase 
(Maraston 2005, Chariot & Bruzual 2007) produce shifts of 
up to ~ 0.2 dex towards lower A4. Finally, the uncertainty 
in the absolute value of the M related to the assumptions 
about the initial mass function (IMF) is within a factor of 
2 for the typical IMFs usually adopted in the literature. 

In this paper, we adopt stellar masses derived with 
Hyperzmass (see Pozzetti et al. 2007), a modified version 
of the photometric redshift code Hyperz (Bolzonella et 
al. 2000). We used smooth exponentially decreasing SFHs 
{SFR{t) oc exp(— t/r) with timescale r = [0.1, oo] and age 
t = [0.1,20] Gyr constrained to be shorter than a Hubble 
time at each redshift). This parametrisation may be inap- 
propriate for actively star-forming galaxies (especially at 
high redshifts, Renzini 2009), because it assumes that all 
galaxies are observed while having their minimum SFR. 
However, for the redshift range explored (z < 1) and for 
passively evolving galaxies this assumption does not ap- 
preciably affect the mass determinations. We checked, for 
example, that the inclusion of delayed exponential SFHs 
{SFR{t) oc t/r^ X exp(— i/r) do not significantly affect the 
stellar masses estimates. As we noted previously, small sec- 
ondary bursts may cause a systematic increase in the mass 
estimate (Fontana et al. 2004). For this study, we adopted 
the Calzetti et al. (2000) extinction law, solar metallicity 
(see also Table 1 in Pozzetti et al. 2007), Bruzual & Chariot 
(2003) population synthesis models, and a Chabrier IMF 
(Chabrier 2003) with lower and upper cutoffs of 0.1 and 
100 Mq, to which we collectively refer hereafter as the 
default parameters, if not specified otherwise. All GSMFs 
were computed by using in addition M05 and CB07 popula- 
tion synthesis models. In particular, we are aware that the 
estimate of ETC stellar masses at high-z are particularly 



sensitive to the TP-AGB phase of the stellar population 
around 1-2 Gyr, producing a shift to lower masses of up to 

0. 2 dex (Maraston et al. 2006). In the redshift range sam- 
pled here (0.1 < z < 1.0), ETGs have on average relatively 
old stellar populations (Thomas et al. 2005) that are not 
dominated by the TP-AGB phase, in particular for massive 
objects. At higher redshifts (z > 1 — 1.5), the effect of the 
TP-AGB phase in younger (1-2 Gyrs) ETGs becomes very 
important and should not be neglected. In our data (lim- 
ited to z < 1), we find that by using M05 models instead 
our main results and conclusions remain almost unchanged 
and we discuss the differences later in the text. 

4. Galaxy classification 

We used 5 different methods to classify each galaxy as ei- 
ther early- or late-type (ETGs and LTGs, hereafter), ac- 
cording to its colour, star formation activity, morphology, 
or spectroscopy, as follows; 

1. Photometric classification (red and blue galaxies): 

We derived galaxy photometric types (VT) from the 
best-fit to the multi-band photometry (from u* to Ks 
band, see previous section). Following Zucca et al. 
(2006), we used the empirical set of (62) SEDs described 
in Ilbert et al. (2006). These SEDs were derived by in- 
terpolating between the four local observed spectra of 
Coleman et al. (1980, ranging from that of the old stel- 
lar population of both M31 and M81 to Sbc, Scd, and 
Im SEDs) and two starburst SEDs from Kinney et al. 
(1996). These templates were also linearly extrapolated 
to the ultraviolet (A < 2000 A) and near-infrared wave- 
lengths using the GISSEL synthetic models (Bruzual & 
Chariot 2003). Galaxies were divided into four types, ac- 
cording to their spectral energy distribution between the 
UV and near-IR. These types correspond to a red E/Sa 
template [VT—l), an early spiral template {VT—2), 
a late spiral template {'PT=3), and an irregular or a 
starburst template (7'T=4). Using the data from the 
VVDS spectroscopic survey, Zucca et al. (2006) verified 
statistically the consistency between the photometric 
classification and the average galaxy spectral proper- 
ties. Using stacked spectra of all galaxies of each of the 
four types, they detected an increasingly bluer contin- 
uum with stronger emission lines from VT=l to 7'T=4, 
confirming the robustness of this classification scheme. 
In Fig. [T] (left panel), we present the U — B rest-frame 
colour distributions of galaxies in our sample divided 
according to the photometric types. From this figure, 
it is evident that the colour bimodality present in the 
sample can be described well by dividing the sample 
according to their photometric type: red galaxies that 
have VT^l (2 103 objects, also called SED-ETGs here- 
after) and blue galaxies that have VT— 2, 3, 4 (6 347 
objects, SED- LTGs), respectively. In Fig. [1] (left panel), 
we also show the colour-mass diagram, which confirms 
that when our photometric classification is applied to 
the lOk-bright zCOSMOS sample it is consistent with 
previous selections of red-sequence galaxies based on 
colour-luminosity diagrams (Bell et al. 2004 and Cimatti 
et al. 2006) and has the advantage of using the entire 
multiband photometric coverage rather than only two 
bands. In principle, given the good multiband coverage, 
the SED fitting is able to break the degeneracy between 
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Fig. 1. Left figure: The upper panel shows the U — B 
rest-frame colour distribution of galaxies with different 
photometric types. The bottom panel shows the colour- 
mass diagram. Right figure: SSFR of galaxies with differ- 
ent photometric types. The upper panel shows the com- 
parison between estimates from the SED fitting and from 
[Oil] for galaxies in the range 0.5 < z < 0.9. The mid- 
dle panel shows the SSFR distribution of galaxies, along 
with the two dotted lines used to separate active from qui- 
escient (log(S'S'Fi?/Gyr^^)= —1) or from passive galax- 
ies {\og{SSFR/GyT-'^)= -2). The bottom panel shows 
the log(S'5'i^i?/Gyr~-'^)-mass diagram. Red dark and light 
cyan colours represent SED-ETGs (VT^l) and SED-LTGs 
{VT=2, 3, 4), respectively. 



a dust-extinguished star-forming galaxy with a red but 
smooth SED and a passive, old galaxy with a strong 
4000 A break in its SED (Pozzetti & Mannucci 2000). 

2. Star formation activity classification (active, quiescent, 
and passive galaxies): 

We divided the sample according to the galaxy star for- 
mation rate {SFR) activity. We used the SFR, and also 
the specific star formation rate {SSFR=SFR/M) de- 
rived from the SED fitting (see previous section) find- 
ing, in general, a good correlation between these and 
the SSFR derived from [Oil] (Maier et al. 2009) for 
galaxies in the redshift range 0.5 < z < 0.9 (shown 
in Fig. [11 right upper panel). We therefore divided the 
sample into active and quiescent galaxies depending on 
whether log(S'S'Fi?/Gyr^^) is above or below —1, i.e., 
galaxies that would take less or more than 10 Gyr to 
double their M at the present SFR, respectively. In our 
sample, we find 5051 active galaxies and 3403 quiescent 
galaxies. With this definition, we note that the num- 
ber of galaxies in the high SSFR sample is ^ 80% of 
the blue galaxies defined on the basis of the SED fitting 
(see Fig. [2 right middle panel). In addition, we investi- 
gate in this paper the population of "passive" galaxies, 
i.e., 1612 galaxies with \og{SSFR/GyT-^)< -2. This di- 
viding value separates well the two parts of the bimodal 
distribution of the SSFR (Fig.[Tl right middle panel). 

3. Morphological classification (spheroidal and disc + ir- 
regular galaxies): 

Taking advantage of the COSMOS HST+ACS images 
(Koekemoer et al. 2007) , we also divided the sample ac- 
cording to morphological class. To take into account un- 
certainties in the morphological classification, we used 
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Fig. 2. U — B rest- frame colours of galaxies with different 
morphological classes using ZEST (left panel) and MRS 
(right panel) classifications. The colour-mass diagrams are 
shown in the bottom panels: red dark and light cyan 
colours represent morpho-ETGs (spheroids) and morpho- 
LTGs (disc-|- irregulars) , respectively. 



two of the available estimates, i.e., the Zurich Estimator 
of Structural Types (ZEST; Scarlata et al. 2007a) and 
a non-parametric estimate derived by our collaborators 
in Marseille (MRS hereafter, Cassata et al. 2007, 2008, 
Tasca et al. 2009). Only a small fraction (~ 3-4%) of 
objects do not have a morphological classification. 
ZEST quantitatively describes the galaxy structure by 
performing a principal component analysis (PGA) in 
the five-dimensional parameter space of asymmetry (A) , 
concentration (G), Gini coefficient (G), the second- 
order moment of the brightest 20% galaxy pixels (M20; 
Abraham et al. 2003), and the ellipticity of the light 
distribution (e). The PGA indicates that the first three 
PG variables account for more than 90% of the variance 
in the original data set, thus almost completely describ- 
ing the galaxy structure. The ZEST classification asso- 
ciates each PG value with a type (ZEST type, ZT, =1 
for early-type galaxies; =2 for disc galaxies; and =3 for 
irregular galaxies) and a "bulgeness" parameter accord- 
ing to the median value of the distribution of Sersic in- 
dices n (Sargent et al. 2007) of all galaxies brighter than 
/ = 22.5 (bulgeness=0,l,2,3, from n > 2.5 being bulge- 
dominated to n < 0.75 being disc-dominated galaxies). 
In this study, the total sample was divided into morpho- 
ETGs and morpho-LTGs, with the morpho-ETGs sam- 
ple including 1680 galaxies classified by ZEST as ellipti- 
cal {ZT—1, 759) or bulge-dominated {ZT=2.0) galax- 
ies, while the morpho-LTGs subsample includes 6413 
galaxies (with ZT greater than 2.0, i.e., from not-bulge- 
dominated to disc-dominated and irregular galaxies). In 
addition, we distinguish ellipticals (ZT—l) from bulge- 
dominated galaxies {ZT— 2.0). 

For comparison, we also used MRS morphological esti- 
mates. This classification scheme separates galaxies on 
the basis of their position in the multi-dimensional pa- 
rameter space of the four non-parametric diagnostics of 
galaxy structures (G, A, G, and axial ratio). First, the 
structural parameters are measured for all galaxies in 
the sample. Then, a randomly extracted subset of 500 
galaxies is visually classified (LT, PG) as ellipticals, spi- 
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rals, or irregulars. This reference catalogue is used to ex- 
plore the multi-dimensional parameter space. For each 
new galaxy that needs to be classified, the distance in 
the parameter space to the 500 reference galaxies is mea- 
sured, and the 11 closest selected. Finally, each galaxy 
is assigned to the most frequent visual class among 
these 11 nearest reference galaxies. This procedure al- 
lows us to convert structural parameters into morpho- 
logical classes: spheroidals (MRS type, MT=1), spi- 
rals {AiT=2), and irregulars (A1T=3). For this paper, 
we also therefore used this morphological classification 
to divide the sample into morpho-ETGs and morpho- 
LTGs, as in the case of ZEST morphologies, yield- 
ing 2383 spheroids {MT—l) and 5831 disc-|-irregular 
{MT=2+3) galaxies, respectively. 
In Fig. m we show the U — B rest-frame colour dis- 
tributions and colour-mass diagrams of these different 
morphological classifications. For both morphological 
classifications, there is a non-negligible number of blue 
spheroidal-type galaxies (~ 37% in MRS and 29% 
in ZEST) (see Sect. 18. 2[) as well as red disc-|-irregular 
galaxies (- 13% in ZEST and - 9% in MRS). 

4. Spectroscopic classification (emission- and absorption- 
line galaxies): 

Using measurements of spectral features (Lamareille 
et al. in preparation), we divided the sample into 
spectral-ETGs and spectral-LTGs following the crite- 
ria defined in Mignoh et al. (2009) for the zCOS- 
MOS spectroscopic sample in the EWo[OII] - D4000 
plane. In the redshift range 0.55 < z < 1.0, 
we find 1079 galaxies without strong emission lines 
and large 4000 A breaks {EWo{[OII]) < 5 A and 
1)4000 + 0.33log{EWoi[OII])) > 1.5, and spectro- 
scopic type ST=1, i.e., spectral-ETGs) and 3304 galax- 
ies with strong emission lines and small 4000 A 
breaks {EWo{[OII]) > 5 A and 1.50 < D4000 -I- 
0.331og(i;Wo([0//])) < 2.22, spectroscopic type 
ST— 2, i.e., spectral-LTGs), respectively. 

5. Combined classification: a clean sample of "bona-fide 
ETGs": 

Given the remaining controversy about morphology 
or colour-selected early-type galaxies (see Franzetti et 
al. 2007 for a discussion of the contamination colour- 
selected samples), we defined a more conservative sam- 
ple of "bona-fide ETGs" by combining different crite- 
ria related to colours, morphology, and spectral proper- 
ties. From the sample of red galaxies (with photomet- 
ric type=l), we selected only those whose properties 
are more accurately described by the 4 reddest tem- 
plates (^ 70%) and we removed galaxies with either 
strong emission hues {EWa{[OII]) or EWo{Ha) > 5 
A, ~ 30%) or rehably classified disc-|-irregular mor- 
phologies (using the intersection between MRS and 
ZEST estimates, ~ 20%) or strong emission at 24 /im 
{K — m{2AiJ,m) > —0.5, ^ 5%), obtaining a final sample 
of 981 red passive spheroids ("bona-fide ETGs"). For 
further details, we refer to Moresco et al. (2010). 

5. The galaxy stellar mass function 

5.1. The method 

To derive the galaxy stellar mass function (GSMF), we fol- 
low traditional techniques used in computing the luminos- 
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Fig. 3. Stellar mass as a function of redshift (small grey 
dots). Also shown are A4\ini (intermediate size dark dots) 
and Mmin (big black circles and lines) at 95% of the M/L 
completeness level (see text): upper panel for the total pop- 
ulation, middle and lower panel for SED-ETGs and SED- 
LTGs, respectively, defined by their photometric type. Also 
shown with a dashed black line is the Ai threshold (cor- 
responding to 95% of the completeness level for the global 
GSMF) derived using mock Millenium samples by Meneux 
et al. (2009). 



ity function. Here we apply the classical non-parametric 
1/V^nax formalism (Schmidt 1968) and estimate the best- 
fit Schechter (1976) parameters {a,Ai*,4>*). In performing 
the fit, we allow ourselves to use up to 2 Schechter function 
components. 

To correct for both the non-targeted sources in spec- 
troscopy and those for which the spectroscopic redshift 
measurement failed, we use a statistical weight associated 
with each galaxy with a secure redshift measurement. This 
weight is the inverse of the product of the target sam- 
pling rate (TSR) and the spectroscopic success rate (SSR). 
Accurate weights were derived by Bolzonella et al. (2009, 
see also Zucca et al. 2009) for all objects with secure 
spectroscopic redshifts, taking into account the magnitude, 
colour, and redshift dependence of the SSR, as well as the 
objects observed as compulsory (^ 2%) and secondary ob- 
jects (~ 2%) in the slits, which have different TSRs from 
the whole sample. 
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5.2. The limits in mass 

In a magnitude-limited sample, the minimum stellar mass 
for which observations were completed depends on both 
the redshift and the stellar mass-to-light ratio Ai/L. This 
latter quantity obviously depends on the stellar populations 
and therefore on galaxy colours. To account for this limit, 
we define at each redshift a minimum mass, A^min, above 
which the derived GSMF is essentially complete because 
all types of galaxies are potentially observable above this 
mass. 

To derive A^min, we calculate the limiting stellar mass 
(A^iim) of each galaxy, i.e., the mass it would have, at its 
spectroscopic redshift, if its apparent magnitude were equal 
to the limiting magnitude of the survey (/um — 22.5) given 
by log(A^iim) = log(A^) + 0.4(7 - h,^). The result is a 
distribution of limiting stellar masses, A^iim, that reflects 
the distribution of stellar A4/L ratios at each redshift in 
our sample. To derive a representative limit for our sample, 
we use the Mum of the 20% faintest galaxies at each red- 
shift. This choice takes into account the colour-luminosity 
relation and therefore includes only galaxies with a typi- 
cal M./L close to the magnitude limit. By doing this, we 
avoid the artificial use of a too stringent limit related to 
the brightest and reddest (with the highest M./L) galaxies, 
which do not significantly contribute close to the magnitude 
limit of the survey. Figure [3] shows the distribution of stel- 
lar masses for all galaxies and of A^um for the 20% faintest 
galaxies. We then define Mrain{z) as the upper envelope 
of the Aliim distribution below which lie 95% of the A^um 
values at each redshift. This A^min corresponds to a 95% 
completeness limit to the M./L ratio at each redshift ob- 
servable by the survey, and is taken to be the completeness 
limit of the GSMF. Meneux et al. (2009) instead used mock 
survey samples to derive a mass limit. For a given redshift 
range and mass threshold, the completeness was simply de- 
fined as the ratio of the number of galaxies brighter than 
the observed flux limit to those of all fluxes. Interestingly, 
even if this method is mo del- dependent, it leads to a similar 
completeness limit as the previous one. For example, for the 
global population we find that our A^min are in fairly good 
agreement with the 95% completeness limit (black curve in 
Fig. El) derived by Meneux et al. (2009) for zCOSMOS, in 
most of tlie redshift range, and in any case the completeness 
for A^>A^i„in is never lower than 85% at any redshift. 

We derived A^min at each redshift and for each galaxy 
subsample used. For the global population, we note that 
A4inin at low redshift is close to the limit of the bluest 
population (see Fig. ^ , because blue galaxies dominate the 
zCOSMOS population close to the magnitude limit at that 
redshift, while the global A^min shifts towards that of the 
reddest population at high redshift (at which redshift a 
significant number of red galaxies is present at the limit of 
our survey, and have a higher Ai/L and therefore higher 
mass limit). 

The GSMFs derived using the 1/Vmax technique, which 
corrects in addition for volume incompleteness, are formally 
complete for A4>A^inin(2inf), i-e., A4min at the lowest red- 
shift of the considered bin. In the parametric fit, we esti- 
mate the best-fit Schechter parameters using data above 
A4„iin(^inf) and plotting as lower limits the data below 
Almin, where the GSMFs are incomplete. In addition, we 
took into account the upper limits that are above the max- 
imum mass found in each considered redshift bin, deriving 



1(7 upper limits by following Gehrels (1986) in the case of 
Poisson statistic for n ~ events (i.e., < 1.84 at Icr). 

6. The bimodality in the zCOSMOS global galaxy 
stellar mass function 

The resulting GSMFs of the zCOSMOS galaxy sample were 
derived for the redshift range 0.1 < z < 1. We divided 
this redshift range into 4 redshift bins to obtain an ap- 
proximately comparable number of objects in each redshift 
bin. This choice also dilutes the effect of the most promi- 
nent large-scale structures at all redshifts (Lilly et al. 2009; 
Kovac et al. 2010). Using the deviations from the median 
number densities in numerous redshift bins, affected by dif- 
ferent structures, we estimate that the systematic error in 
our GSMFs caused by cosmic variance is about 15-20%. 

Figure S] shows the global GSMFs in different redshift 
bins derived using the 1/Knax technique, with the associ- 
ated Poisson errors. The most clearly evident results are 
that at z < 0.55 the global GSMFs, derived using the 
1/^nax technique, are bimodal, exhibit an upturn with a 
steep slope below Al~ 10^-^ A^o, and are not reproduced 
well by a single Schechter function. At higher redshifts, the 
narrower mass range due to the increase in A^min does not 
allow us to explore the upturn in the low-to-intermediate 
mass regime. We therefore allowed the fit to include up to 
two Schechter functions. We verified that the characteristic 
mass (A^*) of the second Schechter function, which domi- 
nates at low mass, is poorly constrained and consistent with 
the characteristic mass of the Schechter function dominat- 
ing at high mass. Therefore, following Baldry et al. (2008), 
we considered hereafter a double Schechter function with a 
single M* value, given by 

«^„a^,..-/-[«(^)"^«(^)"■] 

where (j);^ dA^ is the number density of galaxies with mass 
between A^ and A^ -f dA^, and a2 < ai, so that the second 
term dominates at the lowest masses. 

Using galaxies with A4>Alniin, we find that the slope 
parameters of the 2 Schechter functions differ significantly. 
Using the F-test, we find that the fit with two Schechter 
functions is tighter than the fit with a single function at 
> 3cr confidence level in the first two redshift bins. We 
checked that our findings for the double Schechter shape 
are insensitive to the particular SFHs and set of population 
synthesis models used to estimate A4. This is shown in Fig. 
[5l where we plot the GSMF derived in the first redshift 
bin (0.1 < z < 0.35) using different SFHs and population 
synthesis models to obtain the mass estimates. 

The strength of the bimodality in the global GSMF is 
related to the mutual ratio of the two Schechter functions 
that dominate the fit at the high and low-mass end, respec- 
tively. In Sect. [HI we explore whether the populations of 
ETGs and LTGs can explain this bimodality in the GSMF. 

For many years, there has been evidence of a 'faint-end 
upturn' in the optical galaxy luminosity functions, whereby 
the luminosity function is found to rise steeply about 3~ 
5 magnitudes below the characteristic luminosity (L*), in 
both clusters (Driver et al. 1994, Popesso et al. 2006) and 
the field (Zucca et al. 1997; Blanton et al. 2005). We note, 
however, that the upturn is not evident for all samples 
(Norberg et al. 2002). The existence of a similar upturn 
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Fig. 4. Global galaxy stellar mass functions in four different redshift bins. Red filled points represent 1/Kiiax determi- 
nation and associated Poisson errors, while the lines represent the Schechter fits. The global GSMFs are not reproduced 
well by a single Schechter function (dotted lines), but by two Schechter functions (continuous lines) up to z = 0.55. In 
the first two redshift bins, the GSMFs show the inflection points of the bimodality around Ai^ 5 x 10^ -Mq. The dashed 
vertical lines represent the mass limit in the corresponding redshift bin (A^min)- Data are plotted as lower limits below 
A^min- Upper limits at la (i.e. 1.84 objects) are shown at the high-mass end. Small black dots in all panels represent 
the local GSMF by Baldry et al. (2008). The cross at the right top of the first panel shows an estimate of the cosmic 
variance and mass uncertainties. 



is less clearly evident in the GSMF. Baldry et al. (2004, 
2006, 2008) fitted a similar shape to the local GSMF in the 
SDSS with an upturn at masses below 10^ -Mq- Baldry 
et al. (2008) converted the galaxy luminosity function of 
rich galaxy clusters in the SDSS (Popesso et al. 2006) to an 
equivalent mass function and found an even more promi- 
nent upturn in their GSMFs than that found for the field 
GSMF. This upturn has not been clearly investigated un- 
til now in the GSMFs derived from deep surveys, even if 
a substantial population of low-mass galaxies (< IO^Mq) 
at low redshift {z ~ 0.2), consisting of faint blue galax- 
ies, has been noticed in deep surveys such as the VVDS 
(Pozzetti et al. 2007). We note, a posteriori, by compar- 
ing with other GSMFs from deep surveys (VVDS, Deep2, 
MUSIC, GOODS, PDF, C0MB017) that previous surveys 
(in particular MUSIC, Fontana et al. 2006, and C0MB017, 
Borch et al. 2006) show some hint of bimodality in the shape 
of their GSMFs, even up to z ~ 1 (see Fig. 10 in Pozzetti 
et al. 2007). Ilbert et al. (2010) also found that the GSMFs 
derived in COSMOS using photometric redshifts cannot be 
reproduced well with a single Schechter function, but the 
upturn is less evident than in the zCOSMOS sample. This 
may be caused by the higher precision of the A4 determi- 
nation achieved by using spectroscopic redshifts instead of 
photometric ones. 

In the zCOSMOS survey, the dependence of the GSMF 
shape, on the galaxy environment, as seen in the local 
Universe (Baldry et al. 2006), is investigated in a paral- 
lel paper (Bolzonella et al. 2009), in which the bimodality 
in the global mass function is found to be even stronger in 
high density enviroments. 



In the following sections, we study the evolution of the 
global GSMF and how its shape is related to the bimodal- 
ity in the galaxy properties. We derive the contribution of 
the different galaxy types to the GSMF and its bimodal 
shape, and we explore its evolution with cosmic time. In 
Sect. 17. 11 we also explore how the evolution related to SFHs, 
which differ between high- and low-mass galaxy popula- 
tions, changes the shape of the GSMF bimodality. 



7. The evolution of the galaxy stellar mass function 

We explore how the evolution of the global GSMF with 
cosmic time depends on galaxy mass. From Fig. 31 a mass- 
dependent evolution of the global GSMFs is clearly evi- 
dent: the massive tail is almost stable up to z = 1, while 
the number density of less massive galaxies increases con- 
tinuously with cosmic time, in agreement with previous 
studies (Fontana et al. 2006, Pozzetti et al. 2007). This 
mass-dependent evolution suggests that most of the mas- 
sive galaxies assembled their mass earlier than lower mass 
galaxies ("mass-assembly downsizing"). Figure [6] shows the 
number densities (pw), derived by the 1/V^nax method, for 
two different mass limits of log{M / Mq)> 9.77, 10.77 (cor- 
responding to \og{M/MQ)> 10,11 for a Salpeter (1955) 
IMF, as often used in the literature). We compare our 
results with literature data from other high-z deep sur- 
veys (VVDS, DEEP2, GOODS-MUSIC, FORS Deep Field, 
C0MB017; see Pozzetti et al. 2007 for details and refer- 
ences), and for the local Universe (Cole et al. 2001, Bell et 
al. 2003, Baldry et al. 2008). The local values for zCOSMOS 
were derived assuming an evolution from the lowest redshift 
bin to z = (see next section for details). zCOSMOS data 
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Fig. 5. Global GSMF in the first redshift bin using different 
SFHs and population synthesis models to derive the stel- 
lar masses (BC03, M05, CB07, CB07b=CB07+ secondary 
burst). Points and lines have the same meaning as in Fig,|31 
The dot-dashed lines represent the 2 Schechter components: 
they intersect around 10^° A^o- 



are consistent with most previous observations at 2; < 1 and 
are more accurately determined. We confirm a continuous 
evolution in number density for log(A^/A4Q)> 9,77, which 
increases towards that of the local Universe, and a slightly 
milder evolution for log(A^/A^0)> 10,77, which is negligi- 
ble for z < 0,7 (< 20%) and more rapid above this redshift 
(< 60% since z - 1), 

The same trend with mass is evident in Fig, [3 which 
shows the number densities (pn) as a function of red- 
shift using a more extended range in mass limits (from 
log(A^/7W0)> 9 to \og{M/MQ)> 11,3), Over the en- 
tire mass range, our data are consistent with a faster and 
steeper increase in the number densities with cosmic time 
going from high to low-mass galaxies: the evolution from 
z = 0.44 to z = is 0,12 ± 0.02 dex for log(A^/A^0)> 9.5 
and 0.09 ±0.03 for log{M/MQ)> 10.6. Only galaxies with 
log(A^/Al0)> 11 show negligible evolution from z = 1 to 
the present time (0.03 ± 0.07 dex), while the evolution is 
0.19 ±0.03 dex for log{M/MQ)> 10.6, the lowest mass for 
which we are complete over the whole redshift range. 

zCOSMOS data are consistent with a mass-dependent 
assembly history, with more massive galaxies evolving ear- 
lier than lower mass galaxies (mass-assembly downsizing). 
We therefore confirm that mass is an important parameter 
driving galaxy evolution. In the context of the hierarchi- 
cal scenario, galaxy formation is predicted to be controlled 
primarily by their dark matter halo mass, and therefore, 
one would expect the evolution of individual galaxies to be 
affected also by their enviroment. The relative importance 
of stellar mass and galaxy environment to the shape and 
evolution of the GSMF is investigated for the zCOSMOS 
dataset in a parallel paper (Bolzonella et al. 2009). This 



Fig. 6. Cosmological evolution of the galaxy number den- 
sity as a function of redshift, derived using l/Knax from 
the zCOSMOS for 2 different mass thresholds (> 10^-'^'^ Mq 
and > 10^°'^^ Af© from top to bottom). Dotted errors in- 
clude cosmic variance estimates. zCOSMOS data at z = 
have been derived from the "SFH-evolved GSMFs" (see 
Sect. \7\ and Fig. [H]). The dashed lines correspond to the 
no-evolution solution normalized at z = 0. Results from 
previous surveys (small black points) are also shown. 
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Fig. 7. Number density evolution using different mass lim- 
its for stellar masses. zCOSMOS data at z = have been 
derived from the "SFH-evolved GSMFs" (see Sect. [7] and 
Fig. [8]). The dotted lines correspond to the no-evolution 
solution normalized at z = 0. 
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Fig. 8. Upper panels: "SFH-evolved GSMF^' (empty red points) compared to observed GSMF at redshift z (black filled 
points, from Fig. 2]). In the first panel, we plot the GSMF evolved from the first redshift bin to z = (empty red points), 
compared with the local GSMF (Baldry et al. 2008, shown as small dots). Vertical lines represent the respective mass 
limits. Lower panels: Ratio of "SFH-evolved GSMF" to observed GSMF (both estimated using l/T^nax technique). The 
dotted lines correspond to ±30% around 1. 



study found an accelerated trend of downsizing in over- 
dense regions (see also lovino et al. 2010). 

A detailed comparison of the global GSMF with hi- 
erachical semi- analytical models (SAM) is postponed to 
a future paper. In Sect. O we compare the predictions 
of SAMs with observational ETG GSMFs. Here we recall 
that most hierarchical galaxy assembly models are unable 
to fully account for the observed GSMF and its evolu- 
tion (see Fontana et al. 2004, 2006; Caputi et al. 2006; 
Kitzbichler & White 2007; Marchesini et al. 2009, Fontanot 
et al. 2009 for a detailed comparison with models) . A revi- 
sion of the physical treatment of the baryonic component, 
such as its star-formation history/timescale, of the role of 
feedback, dust content, and/or AGN feedback (Menci et 
al. 2006; Monaco et al. 2007; Bower et al. 2006), and the 
introduction of the uncertainties in the mass determina- 
tion, may help to reduce the disagreement (Kitzbichler &; 
White 2007; Cattaneo et al. 2008; Fontanot et al. 2009). 
Fontanot et al. (2009) compare a broad compilation of avail- 
able data sets with the predictions of three different semi- 
analytic models of galaxy formation within the ACDAI 
framework. When observational errors on stellar mass are 
taken into account, they also find that the models accept- 
ably reproduce the observede number density of massive 
galaxies {\og{M/M(7,)> 11), but that low-mass galaxies 
(log(A^/A^0)= 9 — 10) are predicted to form too early in 
the models and are too passive at late times. Thus, the 
models do not correctly reproduce the downsizing trend 
in stellar mass. We show in Sect. IH] that the typical un- 



certainties in the mass determination cannot fully account 
for the observed excess of massive galaxies relative to some 
hierarchical model predictions, and, when applied, the hier- 
archical models tend to overestimate the high-mass galaxy 
number density at low redshift (see also Fig. 1 in Fontanot 
et al. 2009). 

In the following section, we analyse whether the evo- 
lution of the observed GSMF with cosmic time is driven 
mainly by the SFHs at any given mass, and quantify the 
residual importance of merging events. 

7.1. The importance of the SFHs to the evolution at z < 1 

We attempt to interpret the mass-dependent evolution with 
cosmic time found for the GSMF, i.e., the increase in num- 
ber density at a given M or the shift in at a given 
number density, in terms of only a "pure growth in stel- 
lar mass" with cosmic time related in turn to the mass- 
dependent star-formation history of each galaxy (Thomas 
et al. 2005, Noeske et al. 2007a, b). In this way, we aim 
to quantify the residual importance of merging events or 
other processes after accounting for the role of the galaxy 
star formation activity, which drives the growth in stellar 
mass with cosmic time. A "SFH driven evolution" was also 
adopted by Bell et al. (2007), assuming a constant specific 
SFR between two adjacent redshift bins and an instanta- 
neous return fraction of 45%. 

For each galaxy, we derived the Ai evolved following 
its SFH (exponentially decreasing) derived from the SED 
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fitting. Indeed, the best-fit parameters of the stellar popu- 
lation models (such as the star-formation history timescale 
r and the age) derived from the observed SEDs, allow us to 
follow the evolution of stellar mass in cosmic time, based on 
the assumption that the galaxies evolved in isolation (i.e., 
without merging) and continued to form stars in a smooth 
way with the same SFH derived from the SED fitting. For 
each galaxy at redshift z, we derived the 7W evolved to the 
mean redshift of the previous (i.e., lower) redshift bin. The 
cosmic time elapsed between the mean redshifts of two ad- 
jacent redshift bins ranges between 1.1 and 1.6 Gyr. For 
galaxies in the lowest redshift bin (zmoan — 0.25), we pre- 
dicted the evolved to z = 0, i.e., after about 3 Gyr. We 
performed simulations to test the reliability of this "evolved 
stellar mass" (A^cvoivcd): using models, we created a sim- 
ulated multi-band catalogue at the depth of our catalogue 
(following Bolzonella et al. 2000), from which we were able 
to recover not only the present A4 of the input model (see 
also Pozzetti et al. 2007) but also the model "evolved stellar 
mass" after a given time; the agreement overall is satisfac- 
tory, with no systematic shifts and with the typical disper- 
sions caused by statistical uncertainties and degeneracies 
being of the order of 0.13 — 0.21 dex, depending on the 
elapsed time. 

We find that the mass growth derived from the SED 
fitting in our sample is on average about a factor 1.7 (0.24 
dex) for M=10^ '^ Mq and 1.22 (0.09 dex) for M^10^° 
M^Q, but reaches a factor of 30 for the most extreme star- 
forming galaxies with log(A^/A^0)< 10. In Fig. HI we show 
the GSMF predicted at the average redshift of each bin 
using the galaxies in the higher redshift bin and compare 
it with the measured GSMF at the observed redshift. In 
the first panel, we also show the GSMF predicted at z = 
using galaxies in the first redshift bin (0.1 < z < 0.35). 

First of all, we note a modest evolution from z ~ 0.25 
to z = 0, in particular for low-to-intermediate mass galax- 
ies (log(A4/A^0)< 10) with a shift towards higher masses 
being related to the SFHs. The evolved GSMF approaches 
that observed locally (Baldry et al. 2008) with which it is 
rather consistent for \og{M./ Mq)> 9. 

For the zCOSMOS data set, we find reasonable agree- 
ment at each redshift between the "SFH-evolved GSMF" 
and that observed at the same redshift. Close agreement at 
z > 0.3 was also found by Befi et al. (2007) for COMBO- 
17 by assuming a constant SFH. Our result suggests that 
the GSMF evolution is, therefore, driven mostly by smooth 
and decreasing SFHs, because of the progressive exhaustion 
of the gas reservoir, for example from cold gas accretion, 
rather than by merger events or major bursts. Very similar 
results of a large contribution to the mass growth rate by 
the less massive and star-forming galaxies were obtained in- 
dependently by Vergani et al. (2008) in the VVDS. Figure 
[5] (lower panel) shows the ratio of the evolved to the ob- 
served GSMFs in the first three redshift bins for M>M mm 
of evolved galaxies. The agreement is good particularly 
of intermediate-to-low masses (log(A^/7W0)< 10.6), where 
most of the GSMF evolution occurs and could therefore be 
explained mainly by the growth in mass driven by the SFHs 
in intermediate-to-low mass galaxies. At higher masses, 
galaxies have lower SSFRs and therefore the SFHs are un- 
able to explain the evolution observed in the MF, which in 
any case is quite small (see Fig. [7]). In the next section, we 
explore the role of mergers in the evolution of the GSMF. 
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Fig. 9. Rest-frame U — B colour distribution in a mass- 
limited sample (using the most conservative limit, see text) 
in different redshift bins, as observed (solid black histogram 
for the total population) and predicted using the "SFH- 
evolution" from the previous higher redshift bin (red filled 
histogram). The dotted blue histogram in the top panels 
represent the expected z = colour distribution by apply- 
ing the "SFH-evolution" to the galaxies in the first redshift 
bin. Left panels: assuming no evolution in the dust content, 
i.e., the same dust content derived from the SED fitting at 
the observed redshift; right panels: assuming no dust con- 
tent in the evolved galaxies. 



At all masses, we find differences at most of ^ 40%, 
and even lower than 20% for log(A^/A4Q)< 10.3, between 
the evolved and observed GSMFs. These residual differ- 
ences may be related to different processes. Among them, 
one possibility is an additional SFR, i.e., either more pro- 
longed SFHs than inferred from the SED fitting, a different 
SFH functional form (Renzini 2009), or secondary bursts. 
From the SED fitting, we find that for the blue galaxies the 
timescale (t) distribution for the exponential SFHs has a 
median value of 1 Gyr, even if extended to higher values. 
This median timescale is quite short relative to the value 
inferred from the global SFHs (the compilation of Hopkins 
& Beacom 2006 is consistent with r ~ 3.5 Gyr at z < 1). 
On the other hand, the residual differences could also be 
explained by merger events (see next section) or other dy- 
namical processes that are able to induce a growth in stellar 
mass unrelated to the SFH. 

Here, we also checked that the predicted colour dis- 
tributions for the "SFH-evolved galaxies" were consistent 
with the observed colour distributions at different redshifts. 
Figure [9] compares the expected colour distributions to the 
observed ones at different redshifts for two extreme hy- 
potheses about the dust content: (A) (left panel): no evo- 
lution in the dust content, i.e., assuming the same dust 
content at all redshifts as derived from the SED fitting at 
the observed redshift, or (B) (right panel) no dust content 
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Fig. 10. Panel (a): Growth rate in number density after ac- 
counting for star formation evolution. Different thick lines 
and filled colour points represent estimates at different red- 
shifts (dotted, dashed and sohd for z ^ 0.35, 0.55, and 
0.75, respectively). For comparison, thin lines and empty 
points are the major merger rate estimates by de Ravel 
et al. (2010). Panel (b): Merger per Gyr per galaxy as a 
function of mass. Same lines and points as in panel (a). 
Predictions based on cosmological simulations from Stewart 
et al. (2009), for mergers with mass ratio m/A4 > 0.1, are 
plotted as thin lines and black cross for the lowest and high- 
est redshift bin. Panel (c) Merger per Gyr per galaxy as a 
function of redshift. Thick lines and points represent various 
mass ranges (dotted, dashed and solid for log(A^/A^0)~ 
10.55, 10.85, 11.15, respectively). Predictions from Stewart 
et al. (2009) for m/M > 0.1 are plotted as thin lines 
and black cross for log(A1/A^o)~ 10.3 (dotted line) and 
log(A^/A^0)> 11 (solid line). 



in the evolved galaxies. In both panels, we compare the 
colour distribution for galaxies in a mass-limited sample 
using the most conservative limit Almin(-Zsup) for evolved 
masses. Model A predicts too red colours at all redshifts, 
while only the extreme case of model B, i.e., without any 
dust content, is consistent with the observed colour distri- 
butions, at least at z > 0.35, while, in any case, it predicts 
too red colours at lower redshift. This could be related 
(again) to a too rapid decline in decreasing (quenching) 
of the SFRs of each galaxy as inferred by the SED fitting 
technique. We therefore conclude that more extended SFHs 
or secondary bursts appear to be necessary to more accu- 
rately reproduce the observed colours, in particular at low- 
redshift. One possibility would be to directly model the 
observed decrease/quenching of the SFR at later cosmic 
times for progressively lower masses (see also Sect. 01). We 
postpone the discussion of this possibility to a future paper. 



7.2. A limit to the galaxy merging 

The small amount ofevolution observed in the GSMF ar- 
gues against a dominant contribution of galaxy mergers to 
galaxy evolution. If merger events, indeed, were efficient in 
forming the most massive galaxies at z < 1 (as predicted by 
hierarchical models) , we should observe instead a detectable 
increase in the number density of massive galaxies. In par- 
ticular, the negligible evolution in the GSMF for galaxies 
with \og{M/MQ)> 11 up to z ~ 1 suggests that these 
galaxies formed and assembled their stellar mass at higher 
redshifts, ruling out a major role for mergers at z < 1. 

We derived the contribution of mergers to the evolution 
of the GSMF. We first estimated the growth rate of the 
number densities for different M, after accounting for the 
growth in mass produced by ongoing star formation (see 
previous section) . The residual growth rate of the number 
density was calculated to be 



Pn{M,z) 



(2) 



where pN was evaluated using the 1/Knax method. The re- 
sult is shown in Fig. [TU] {panel a) as a function of A4 in the 
various redshift bins. We compare these values with the 
major- merger rate {pmm) derived from galaxy pairs by de 
Ravel et al. (2010) for the zCOSMOS sample. It is evident 
that major-merger rate is always below the growth rate by 
a factor of 2 or more, except in the intermediate redshift 
range (z ^ 0.5). This suggests that major mergers do not 
contribute significantly to the residual evolution. We, there- 
fore, derived the merger rate (major -|-minor) that could ex- 
plain the residual evolution. We estimated the effective in- 
crease in number densities at any given A4 by accounting 
also for galaxies that at the same time leave the current 
mass bin due to major merging. This estimate repre- 
sents the number of galaxies produced by mergers 
(major+minor) per Gyr per Mpc'^ to explain the 
GSMF growth at that mass. Finally, assuming that 
a galaxy is the product of the last single merger 
(within galaxy pairs of any mass ratio m/A4), we 
derived the number of mergers per galaxy per Gyr 
dividing by the number densities of galaxies. We 
note, indeed, that the time interval between two 
adjacent redshift bins is not long (1-1.5 Gyrs) rela- 
tive to the typical merger timescale used in SAMs 
(De Lucia et al. 2010) or predicted by numerical 
simulations (Boylan-Kolchin et al. 2008) that pre- 
vent multiple minor mergers predominating. We, 
therefore, derive the merging rate per galELxy per 
Gyr (/m) in the following way: 



fuiM) 



[pNjM, z) + PmmjM, Zsup)] 

Pn{M,z) 



merger /gal /Gyr. (3) 



The results are shown in Fig. [TU] as a function of mass 
(panel b) and redshift [panel c). At high redshift, we find 
for massive galaxies {\og{AA / Mq)> 10.6 at z > 0.5) that 
the merger rate is quite high, ^ 0.2 — 0.4 merger Gyr"-'^ 
gal~^, but decreases rapidly to ~ 0.1 with cosmic time, be- 
ing only marginally consistent with an increase with mass. 
At lower redshifts, the merger rate does not show any clear 
dependence on mass and always remains below 0.2 at all 
masses. On average, we measure 0.16 merger Gyr~^ gal~^ 
since z = 1. The major merger rate per galaxy derived by 
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Fig. 11. Galaxy stellar mass function by galaxy types. Left panel: GSMFs by photometric types {Vl'=l as red squares, 
PT=2+3+4 as blue triangles, total population as empty black circles). Right panel: GSMFs for morphological types 
(-ZT=1,2.0 as magenta squares and ZT>2.0 as dark green triangles) In both panels, points represent the 1/Vmax deter- 
mination, while continuous lines the Schechter fits. Dotted lines reported in each panel, as a reference, are the Schechter 
fits to the first redshift bin. Dashed vertical lines represent the mass limit for the corresponding redshift bin (A^min)- 



De Ravel et al. (2010) is always below 0.1 at all redshifts 
and masses explored (~ 0.06 on average). Integrated over 
cosmic time (from z ~ 1 to z ~ 0.1), massive galaxies 
{\og{M/MQ)^ 10.6) have experienced about 0.7 merger 
gal~^ (< 0.3 below z ^ 0.7), of which less than 0.2 are 
major mergers. In Fig.[TU]we compared our estimates with 
the model predictions of Stewart et al. (2009) for merg- 
ers within galaxy pairs of mass ratio greater than 0.1 
{m/Ai > 0.1). Galaxy mergers in a hierarchical ACDM 
scenario increase with stellar mass, redshift, and decreasing 
m/M ratio. Minor mergers contribute considerably more 
than major mergers, dominating in all but the most mas- 
sive galaxies (Parry et al. 2009). Indeed, even in hierarchi- 
cal formation models, most of ellipticals and spiral bulges 
acquire their stellar mass through minor mergers or disc 
instabilities. Given our uncertainties, we cannot draw any 
firm conclusions about either the redshift or mass depen- 
dences, even if our data appear to exclude a strong depen- 
dence on mass, at odds with hierarchical predictions. There 
is, instead, some evidence of an increase with redshift, qual- 
itatively in agreement with models. 

We therefore conclude that the role of major merg- 
ing events is not dominant and cannot completely ex- 
plain the evolution of the GSMF, even after accounting for 
the growth in mass produced by ongoing star formation. 



Mergers in general (major -I- minor) at a rate of between 0.1 
and 0.4 Gyr~^ gal~^ may account for the residual evolu- 
tion, these rates possibly decreasing with cosmic time but 
not having a clear dependence on mass. Oesch et al. (2010) 
discuss in greater detail the role of merging in the evolu- 
tion of massive galaxies within the COSMOS survey. In a 
future paper on the final zCOSMOS sample, we will explore 
in detail the role of both major and minor merging on the 
evolution of the GSMF. 



8. Mass functions by galaxy types 

To explore the bimodality observed in the GSMF and how 
it depends on both time and galaxy-type evolution, we de- 
rived a type-dependent GSMF using the different classifi- 
cations described in Sect. 01 

Figure [TT] (left panel) shows the GSMF divided 
into different photometric types, red {VT=1) versus blue 
(■PT=2,3,4) galaxies; the GSMFs obtained by using mor- 
phological criteria to differentiate between ETGs and LTGs 
are shown in the right panel of the same figure. We find that 
the red/ spheroidal and blue/disc-|-irr populations have very 
different GSMFs, which can be linked to the bimodality of 
the global GSMF. The GSMFs of both ETGs and LTGs are 
reproduced quite well by a single Schechter function. The 
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LTG GSMFs exhibit a steep low-mass end, which in con- 
strast with is flat for the ETGs at all redshifts. However, 
we find a hint of a small upturn in both the GSMFs of 
ETGs for the morpho-ETGs (see Fig. [TT|) and in the sec- 
ond redshift bin of our other ETG classifications (see Fig. 

Drory et al. (2009), instead, found that neither their 
red (passive) nor their blue (star-forming) galaxy stellar 
mass functions could be fitted well with a single Schechter 
function, but show an upturn at low masses. In general, 
for all the criteria adopted, we find that the massive end 
of the mass function {Ai> 10^°'^ -M-q) is mainly dom- 
inated by ETGs (red/spheroidal/passive) up to z — 1, 
while LTGs (blue/disc-|-irr/active) mostly contribute to the 
intermediate/low-mass part {M<M*) of the mass function 
at all redshifts. For log(A^/A^0)> 11, the number density 
of blue galaxies is always below ^ 2x 10"'* gal/Mpc'^/logTW. 

Using the 1/Knax data points of the GSMFs of ETGs 
and LTGs, we derived the intersection of the two popula- 
tions (A^cross)- For the photometric classification, A^cross 
is at ~ 1.9(±0.2) X 10^" Mq in the first redshift bin and 
evolves with redshift, increasing by about a factor of 2 
(see red circles in Fig. [T2|) by z = 1. Although for some- 
what higher values of cross, a similar trend with red- 
shift was found using the morphological classifications to 
divide the sample into spheroids and disc-|-irregular galax- 
ies, for both classification schemes explored in this pa- 
per. We determined similar values of cross by dividing 
the sample in terms of the star formation activity (ac- 
tive versus passive with log(S'5'FiZ/Gyr^^)= —2), while 
A( cross is quite low when quiescient galaxies are considered 
{\og{SSFR/Gyi~^)= —1). Our results are consistent with 
previous determinations of cross in other deep surveys, 
such as VVDS (Vergani et al. 2008) and DEEP2 (Bundy et 
al. 2006), the only exception being the value of A^cross de- 
rived by Bundy et al. (2006) using morphological types (the 
highest line in Fig. [T^ . In this case, the difference may be 
caused by a different and more extreme definition of mor- 
phological ETGs by Bundy et al. (2006). The extrapolation 
of our A^ cross to z = is also consistent with local results 
(Baldry et al. 2004; Bell et al. 2003). 

The evolution of A^ cross is indicative of an in- 
crease/decrease in the fraction of early/late galaxies with 
cosmic time at intermediate masses, as already noted in 
previous surveys (Fontana et al. 2004, Bundy et al. 2006, 
Arnouts et al. 2007, Scarlata et al. 2007b, Vergani et al. 
2008). Oesch et al. (2010) found similar results within the 
COSMOS survey using extremely robust photometric red- 
shifts. Bolzonella et al. (2009) found a more rapid evolution 
of A^cross in high density environments (see also lovino et 
al. 2010 and Kovac et al. 2009 for groups). However, previ- 
ous surveys have not clearly established the precise evolu- 
tion in absolute number densities of the two populations, in 
particular for late-type galaxies. By comparing the zCOS- 
MOS GSMFs for different redshift bins, we find that the 
evolution of A^ cross is caused mainly by a clear increase 
with cosmic time in the number density of red/spheroids of 
intermediate mass (around A4* and below 10*^ A^o), cou- 
pled with only a marginal evolution in the number density 
of blue/spiral galaxies at the same masses. 

These trends become even more distinctive when plot- 
ting the number densities, derived using the 1/Vmax 
technique, versus redshift in different mass ranges 
(log(Al/AlQ)= 10 - 11.5), as shown in Fig. M We find 
that the number density of massive (log(A(/A(0)> 11) red 
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Fig. 12. The intersection of the GSMFs (Af cross) of the 
two populations (ETGs and LTGs) as a function of red- 
shift. Values were derived from the l/Knax estimate of the 
GSMFs. Different symbols (and colours) refer to the differ- 
ent classification schemes adopted (see figure legend), up- 
per limits are shown where Atcross is lower than Atniin(z). 
Black crosses show data from the literature (references are 
marked, along with the respective classification criteria). 



types is almost constant up to z = 1 (see also Pozzetti et 
al. 2003, Scarlata et al. 2007b, Ilbert et al. 2010) and dom- 
inates the total number density at these masses, while the 
number density of intermediate-mass (10*° — 10** A^©) red 
types increases with cosmic time. This increase is steeper 
at lower masses. Between z ~ 0.45 and z ~ 0.25, it is 
0.05 ± 0.08, 0.14 ± 0.04, and 0.24 ± 0.05 dex for the three 
considered mass ranges of log(AI/A(0)= 11, 10.5, and 10, 
respectively. An independent study (Brown et al. 2008), 
using halo occupation distributions, shows that while very 
massive halos often double in mass over the past 7 Gyr, 
the stellar masses of their central galaxies typically grow by 
only 30%. Our result is indicative of a mass-assembly down- 
sizing as already noted for the global population and even 
more clearly in the build-up of the red sequence (Cimatti 
et al. 2006). 

In contrast, we find that the number density of late blue 
types with Ai> 10**^ A^© remains approximately constant 
with cosmic time from z = 1. Their evolution in number 
density is less than 10 — 20% over the entire redshift and 
mass range fFig. [Til upper panel). These results are consis- 
tent with those of Arnouts et al. (2007), who found that 
the integrated stellar mass density of the active population 
shows only a modest mass growth rate, in contrast to an in- 
crease by a factor of 2 for the quiescent population. Bundy 
et al. (2006) found that the abundance of blue galaxies de- 
clines by 0.1 - 0.2 dex from z = 0.75 - 1.0 to z = 0.4 - 0.7 
at log(Af/A(©)= 10.6 - 11.3 (see their Fig. 6). A small 
evolution of the "blue cloud GSMF" was also detected by 
Bell et al. (2007) in the COMBO-17 survey out to z ~ 0.9. 

How is this small, possibly negligible evolution in the 
blue GSMFs related to the global decrease with cosmic time 
in the SFR density since z ~ 1 (see Hopkins et al. 2006) 
and the downsizing in the SFR (Cowie et al. 1996)? To 
answer these questions, we analysed the GSMFs of galaxies 
with different SFR activities. In particular, we derived the 
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Fig. 13. Number density evolution in different mass 
ranges. Different symbols and colours refer to different pop- 
ulation classes: black squares refers to tlie total population, 
empty red and filled blue triangles to SED-ETGs and SED- 
LTGs using photometric classification, while filled magenta 
and empty green circles refer to quiescent and active popu- 
lations, respectively. Hatched regions refer to GSMF com- 
plete in the mass range and redshift considered, while dot- 
ted lines and lower limits are plotted at redshifts where the 
GSMF is not complete in the considered mass range. 



GSMFs for the most extreme star-forming population, i.e., 
those with high SSFR. 

Using estimates of the SFR and SSFR derived from 
the [Oil] lines (Moustakas et al. 2006 calibration), Maier 
et al. (2009) find that very few zCOSMOS galaxies with 
masses above log(A^/A^0)> 10.8 have 1/SSFRs below the 
age of the universe at 0.5 < z < 0.7, while there ex- 
ist several dozens of galaxies above the same mass that 
are strongly star-forming at 0.7 < z < 0.9 (see Fig. 2 
in Maier et al. 2009). The same trend is observed using 
the SFR and SSFR derived from the fit to the multi- 
band photometry with synthetic models. Using an SED fit- 
ting determination, we note, first of all, that the median 
SFR and SSFR decreases with cosmic time in massive blue 
galaxies with log(A^/A^0)> 10.5 (for example: the median 
\og{SSFR/Gyi-^) is -1.0 at 0.7 < z < 1.0, while it is -1.4 
at z < 0.3). 

We therefore derived the GSMFs using the SSFR de- 
rived from the SED fitting and by dividing the sample as 
described in Sect. El i.e., as those defining active and qui- 
escent galaxies with \og{SSFR/Gyi~^) above or below — 1. 
FigurefHl flower panel) shows indeed a decrease with cosmic 
time in the GSMFs of the high-SSFR galaxies, in particu- 
lar for log(A^/A40)> 10.3. Their GSMF shows, indeed, a 
mass-dependent evolution that is stronger at the massive 
end. In the highest redshift bin (0.75 < z < 1), the GSMF 
of the high-SSFR galaxies increases in number density and 
approaches that previously derived for the sample of blue 
(SED-LTG) galaxies. 

By comparing the GSMFs of the blue galaxies with 
those of high-SSFR galaxies, we find that, while the num- 
ber density of blue galaxies overall remains approximately 
constant, it decreases significantly with cosmic time for 
high-SSFR massive objects (see Fig. [H)) . From z ^ 0.85 to 




Fig. 14. GSMFs for LTGs using two classification schemes. 
Upper panel: blue galaxies (7'T=1), and lower panel: active 
galaxies {log{SSFR/Gyr^^)> -1), as defined in Sect, g) 
Only 1/Vinax determinations have been plotted, with their 
uncertainties. Different lines and colours refer to different 
redshift range (see legend). 

z ^ 0.25, the number density of high-SSFR active galax- 
ies with log(A^/A^0)> 10.5 decreases by a factor of ~ 3 
(—0.49 ± 0.12 dex). Therefore, the decrease with cosmic 
time in the number of active massive galaxies is balanced 
by the constancy in intermediate-activity blue galaxies and 
the increase of intermediate-mass quiescent and red galax- 
ies. Studying the B-band luminosity function in zCOSMOS, 
Zucca et al. (2009) found results consistent with a scenario 
in which some blue galaxies are transformed into red galax- 
ies with increasing cosmic time. 

To summarize, these data suggest that with cosmic time 
we are witnessing a transformation from active to passive 
galaxies and a corresponding decrease (increase) in the frac- 
tion of late (early) types. These changes in galaxy fractions 
are caused mainly by the clear increase with cosmic time in 
the number density of intermediate-mass (7W~ 10^" — 10^^ 
A4q) early- type galaxies, while the density of intermediate- 
mass blue or morpho-LTGs remains almost constant, be- 
ing continually replenished by blue active galaxies of even 
lower masses. The median SSFR of blue massive galaxies 
decreases with cosmic time. Therefore, blue highly star- 
forming (of high SSFR) galaxies of intermediate mass 
{\og{M/MQ)= 10 — 11) increase in mass but decrease in 
either SFR or SSFR with cosmic time, transforming into 
low- SSFR blue objects and, after the quenching of their 
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SFR, into red passive objects of intermediate mass. Less- 
massive, blue, active objects (log(A^/A^0)<lO) increase in 
mass to replace in the GSMFs the blue intermediate-mass 
objects, whose density consequently remains almost con- 
stant with time. In agreement with our results. Bell et al. 
(2007) found that, if there were a growth in stellar mass 
with a constant SFR equal to the instantaneous one (af- 
ter subtracting 45% of the return fraction), the GSMFs of 
blue galaxies would be dramatically overproduced. By only 
assuming that all the growth in stellar mass is added to 
the red sequence, they reproduced the evolution of the blue 
and red stellar mass functions with remarkable accuracy. 
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Fig. 15. Evolution of the growth rate in number density 
(top panel) as a function of redshift for red galaxies (VT—l) 
in different mass ranges (represented by different symbols 
and colours) . Values have been plotted at the intermediate 
redshift between two adjacent redshift bins. In the bottom 
panel, we show the evolution with redshift of the fraction 
of blue galaxies {VT— 2,3,4) that are transforming into red 
galaxies per Gyr (see text). In all the panels, hatched re- 
gions refer to GSMFs that are complete in the mass range 
and redshift considered, while dotted lines are connecting 
points at redshifts where the GSMF is not complete in the 
considered mass range. 



8.1. The build-up of the red sequence 

We attempt to quantify the growth rate (in terms of both 
number and mass density) of red galaxies as a function 
of redshift in different mass ranges. This rate can be in- 
terpreted as the flux required to describe the migration of 
galaxies with cosmic time from the blue cloud to the red 
sequence. 

Figure [TS] shows the evolution with redshift for the 
red population (VT—l), in different mass ranges, of the 
growth rate in number density {pN, upper panel), es- 
timated using 1/T4iax data points. For log(A^/A^0)< 



11, the growth rate is in the range of 10 ^ - 10 
gal/Mpc^/Gyr/logA^ and is even lower approaching zero 
for higher masses. The growth rate in mass density is 
PM— lO*" - 10^ A^0/Mpc^/Gyr/logA^. Integrating above 
log(A4/X0)=9.8, we find that /9Ar=6.8(±1.2) x 10"'' 
gal/MpcVGyr and pM^ 1.8(±0.4) x lO'' 7W0/MpcVGyr 
between the second and the first redshift bin (z ^ 0.34). 
Despite the large errors, we find a clear increase with cosmic 
time in the growth rate of galaxies with log(A^/A^0)=lO.4- 
10.7, but a hint of a decrease with cosmic time at higher A4 
(log(A4/A40)=lO.7-ll). For log{M/MQ)> 11, the growth 
rate in number density is consistent with zero for almost all 
redshift bins. These trends can be interpreted as a mass- 
assembly downsizing signal, i.e., most massive red galaxies 
assembled their mass earlier than lower mass red galaxies. 
For comparison, Arnouts et al. (2007) for the VVDS found 
that the median mass growth rate (integrated over the en- 
tire mass range log(7W/A^0)=8 — 13) between z ^ 2 and 
z = is 1.7(±0.4) X 10^ Mq/Mpc^/Gyt, which is consistent 
with our values for log(M/MQ)> 9.8. Walcher et al. (2008) 
derived a somewhat higher value of pm= 6.5 — 10 x 10'' 
A^0/Mpc"^/Gyr from z ~ 1 to z 0.5 integrated over the 
whole mass range. The difference from our results, indeed, 
could be due to their more extended mass range. 

In principle, by comparing this growth rate with the 
number density of blue galaxies we should be able to esti- 
mate the fraction of blue galaxies that, at any given time, 
are transforming into red galaxies (/g) after the quench- 
ing of their SFR. We therefore derived the fraction rate of 
quenching 



/q(-M) 



Gyr- 



(4) 



using the 1/Vmax estimation for p^r, as shown in the bot- 
tom panel of Fig. [T5] as a function of redshift and mass. 
For \og{Ai/ A4q)> 11, the rates were unconstrained be- 
cause of the low number density of blue galaxies. The val- 
ues obtained for log(7W/7W0)< 11 have a median value of 
^25% Gyr~^. As for the growth rate, we find that with cos- 
mic time this fraction rate decreases for high-mass galax- 
ies, but increases for intermediate-mass galaxies, at least 
for z < 0.7. At z < 0.4, we are unable to explore the 
transformation from blue to red for the low mass regime 
{\og{M/M(7))< 10.4) and to verify whether the trend of 
the blue fraction increasing with cosmic time for low-mass 
galaxies is shifted to even lower redshift, as expected by the 
downsizing scenario. We note however that all fraction val- 
ues remain very uncertain and their uncertainties increase 
yet further if we include the effects of cosmic variance. 

Vergani et al. (2009) explore the properties and number 
densities of the post-starburst (PSB) galaxy population, 
spectroscopically identified in the zCOSMOS sample, as the 
possible link population in the transition phase between 
the blue cloud and the red-sequence. They found that this 
galaxy population, which is affected by a sudden quenching 
of its star-formation activity, may increase the stellar mass 
density of the red-sequence by up to a non-negligible level 
of 10%. 

8.2. Elliptical/early type galaxy evolution 

Hierarchical semi- analytical models (De Lucia et al. 2006) 
predict that stars in more massive elliptical galaxies, rel- 
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Fig. 16. GSMFs for ETGs using various classi- 
fication schemes: red galaxies {VT—l), passive 
{\og{SSFRlGyY^^)< -2), spheroidals (Zr=l+2.0, 
ZT=1, MT=l), absorption-line galaxies {ST~l), and 
"bona-fide ETGs" as defined in Sect. Sj Only 1/Knax 
determinations have been plotted, with their uncertainties. 
Different lines and colours refer to different redshift ranges 
(see upper left panel). Grey hatched region, bounded by 
solid black lines, refer to local GSMF for ETGs derived 
from Baldry et al. (2008) and Baldry et al. (2006) (see 
text) . 



ative to their lower-mass counterparts, are older but as- 
semble their mass by means of mergers later in cosmic 
time. For example, 50% of ellipticals with \og{A4/MQ)> 11 
form most (80%) of their stars at Zform — 1-6, but as- 
semble their mass by ^assembly— 0.2, which is defined as 
the redshift at which most (80%) of the stars that make 
up the galaxy at redshift zero are already assembled into 
a single object. Galaxies with \og{M/M(?))> 9.6, instead, 
have a median Zform — 1-3 and median Zasscmbiy — 0.9 (see 
Figs. 4 and 5 in De Lucia et al. 2006). As pointed out by 
Cimatti et al. (2006), this trend is the opposite of that of 
the mass -assembly downsizing found for red galaxies in the 
luminosity function, confirmed by Scarlata et al. (2007b), 
i.e., most massive red galaxies assembled their mass earlier 
than lower-mass red ones. In the following, we explore in 
detail for the zCOSMOS data set the evolution of elliptical 
and all early-type galaxies, using the different classification 
schemes described in Sect. |3]based on colours, morphology, 
spectral features, or a combination of them. 



Figure \W\ of the GSMFs for the different popula- 
tions of ETGs: red galaxies {VT=1), passive galax- 
ies (log(S'S'Fi?/Gyr-i)< -2), spheroidals {MT=1, 
ZT=l+2.0, ZT—1), absorption-line galaxies (iS7~=l), and 
"bona-fide ETGs" as defined in Sect. [3 

First of all, we note that there are significant differences 
in the normalization of the GSMFs at all A4, the pure 
ellipticals classified by ZEST (ZT—l) having the lowest 
value at all redshifts and mass ranges explored. In contrast, 
the shapes of the various GSMFs are all similar, there being 
a decline in number density for Ai<Ai* for all GSMFs that 
is described well by a single Schechter function, except for 
the GSMF of the spheroids with A4T—1, which shows an 
excess at low-mass (log(A^/7M0)< 9.5) in the first redshift 
bin. 
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Fig. 17. Stacked spectra for galaxies with MRS elliptical 
morphology (MT—l) and blue colours {VT=2+3+i) at 
low-mass {\og{M/M(7))< 9 upper spectrum in blue) and 
high-mass (log(A^/A^o)> 10.5 lower spectrum in red). 

We then studied in detail this population of low-mass 
spheroids {AiT=l), which dominates the GSMF and ex- 
hibits blue colours {PT >= 2, i.e., {U - B)rest < !)• We 
found 89 of these objects with log(A^/A^0)< 9 at low 
redshift (z < 0.5). They have quite compact sizes, with a 
small median effective radius {r^ — 0.6 kpc). We also iden- 
tified a contribution from blue MRS spheroids {VT>— 2, 
MT= 1) also at high mass {\og{MI Mq)> 10.5) at all 
redshifts (182 objects). This latter population of massive 
blue MRS spheroids does not dominate the high-mass end 
of the GSMF 20%), as the former does the low- mass 
end. Therefore, this suggests that the shape of the GSMF 
for this population of blue spheroids differs from that for 
red spheroids, at least when using MRS morphologies. This 
is consistent with the conclusions of Zucca et al. (2009) 
about the contribution of this population to the i?-band 
luminosity function. However, for these objects, in particu- 
lar at low mass, we note that MRS morphologies do not 
always agree with ZEST ones, which implies that there 
are significant uncertainties in deriving morphologies for 
faint low-mass objects. We, therefore, verified their mor- 
phologies visually (i.e., PC, LT, LP). We found that most 
of the blue low-mass MRS spheroids are bulge-dominated 
spirals (Sa). Most of the massive ones, even if their mor- 
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phological classifications are in agreement between the two 
methods (ZEST and MRS), appear to be spirals, between 
Sa and Sc, although they have morphological parameters 
(asymmetry and concentration) that are more similar to 
those of elliptical galaxies. Composite spectra of this pop- 
ulation were generated (following Mignoli et al. 2009) by 
averaging all their spectra at z < 0.5 and dividing them 
into low-mass and high-mass objects. In Fig. I17[ the av- 
erage spectra are plotted: low-mass blue spheroids show a 
blue continuum and emission lines typical of star-forming 
galaxies, indicating that these galaxies are experiencing sig- 
nificant star-formation, while the high-mass blue spheroids 
have on average a red spectrum with absorption lines that 
are largely indistinguishable from those of the purely pas- 
sive galaxies (see Mignoli et al. 2009), but in addition ex- 
hibit emission lines that could be indicative of some nuclear 
activity contamination (i.e., [NII]/Hq and [OIII]/H^ are 
consistent with LINERs). An AGN component may also 
explain their typically high concentration index obtained 
from their HST images. These AGN might play a role in 
quenching the SFR in this population of blue spheroidal 
massive galaxies, which are possibly in the transition phase 
between the blue cloud and the red-sequence. 

Finally, we studied the evolution with cosmic time of 
the GSMFs of each of the classes of red, passive, ellipticals, 
spheroids, absorption-line, and "bona-fide ETG" galaxies. 
In Fig. [Til we show the GSMFs of these various ETG 
classes at increasing redshift. We also show, as reference, 
the local GSMF of ETGs, which we derived using the 
global GSMF by Baldry et al. (2008) and the fraction of 
red galaxies estimated by Baldry et al. (2006) for their 4 
central environment bins. This is indicated in Fig. [16] by 
the grey shaded region. We note very similar evolutions in 
the ETG population for all the different means of classify- 
ing an ETG. As claimed for red galaxies, we find a clear 
"mass-assembly downsizing" evolution, with no evidence of 
an increase with cosmic time in the number density of the 
most massive ETGs since z ~ 1 and a progressive increase 
in the GSMFs at intermediate mass (log(A^/A^0)< 11) 
out to 2; = 0. Using M05 models to estimate the stellar 
masses, we also find that the same trend is visible in our 
data, and that the use of these models does not strongly 
affect our main conclusions. This is because in the redshift 
range studied here (0.1 < z < 1.0), ETGs have on average 
old stellar populations that are not dominated by TP- AGE 
stars particularly in massive objects. 

For massive ETGs, similar results were first noted in 
the near-IR (if -band) luminosity function by Pozzetti et al. 
(2003) and the mass function by Fontana et al. (2004) in 
the K20 survey, using spectral classification, and later con- 
firmed by larger optical and near-IR surveys, such as the 
VVDS, using colours or spectra to define ETGs (Zucca et 
al. 2007; Arnouts et al. 2007; Vergani et al. 2008). Scarlata 
et al. (2007b) found that in COSMOS both morphologically 
and photometrically selected subsamples of ETGs show no 
evolution in their number density at the bright end of the 
B-band luminosity function {L > 2.5L*) out to 2 ~ 0.7, and 
there is a deficit of a factor of about 2 — 3 of fainter ETGs 
over the same cosmic period, as also confirmed by Zucca 
et al. (2009) in the zCOSMOS sample. Ilbert et al. (2010), 
selecting ETGs (by morphology or colour) in SCOSMOS 
(Sanders et al. 2007) with accurate photometric redshifts 
(Ilbert et al. 2009) up to Zphoto = 2, found similar re- 
sults for massive ETGs at z < 1 (evolution < 0.2 dex for 



log(A^/A^0)> 11) and a more rapid evolution at higher 
redshifts (by a factor of 15-20 between z = 1.5 — 2 and 
z = 0.8 — 1) or lower masses (increasing by a factor of 4.5 
between z = 0.8-1 and z = 0.2-0.4 for log(7W/A^Q)~ 10). 

8.3. The building redshift of ETGs 

Our results suggest that lower mass ETGs (regardless of 
the method used to classify them) assembled their mass 
later than higher-mass ETGs. We can, indeed, estimate the 
redshift at which 50% of ETGs have already built-up their 
mass (zbuiiding), as a function of stellar mass, i.e., at which 
redshift their number density has decreased by a factor of 
2 from z = {(t>z/4'o = 0.5). Since we consider the number 
densities of ETGs only above z > 0.25, we can derive only 
upper limits to ZbuiWing, as defined above, at least in the 
low mass range. For this reason, we used as a local reference 
the GSMF for ETGs, described in the previous section and 
shown in Fig. [IS] as a grey shaded region. 




Fig. 18. Redshift at which the ETG GSMF has decreased 
by a factor 2 {(j)z/<Po — 0.5, i.e., Zbuiiding, see text) as 
a function of M. Left panel: Circles magenta points re- 
fer to data derived within zCOSMOS for passive galaxies 
(log(55'FiZ/Gyr-i)< -2): open circles have been derived 
completely from zCOSMOS dataset, using as a reference 
the GSMF in the first redshift bin (0o = <f>{z 0.24), 
along with the region covered using various ETG classifica- 
tion criteria; filled circles are ZbuiWing using the local ETG 
GSMF (see text), along with the spread in its shape re- 
lated to their enviroment. Dotted magenta curve shows the 
A^min above which we are formally complete. Also plotted 
are zbuiiding from SAMs: WDL08 (shaded region bounded 
by solid blue hues), S08 (shaded region bounded by dot- 
ted green lines), and MORGANA (shaded region bounded 
by dotted black lines), derived using the original GSMF or 
convolved with 0.25 dex in mass uncertainties (see Fontanot 
et al. 2009). Right panel: Zbuiiding for zCOSMOS data points 
are compared with the trend with stellar mass of the me- 
dian Zassombiy and Zform (80% of Star assembled and formed, 
respectively) from De Lucia et al. (2006, DL06) and with 
the redshift at which the stellar formation ends (80% of star 
formed) in local ETGs estimated by Thomas et al (2009, 
T09). 
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Within the zCOSMOS redshift range, we find for ETG 
galaxies that Zbuiiding increases with A4 regardless of clas- 
sification method, with Zbuiiding'^ 0.4, 0.5, and 0.8 for 
logiM/Me)'^ 10, 10.4, and 10.8, respectively (Fig.[lHl). At 
higher mass, the number density of ETGs is almost con- 
stant up to redshift ^ 0.7 and decreases by less than a 
factor of 2 after z ~ 1, thus we set 1 as the lower limit to 
-Zbuiiding for log{M/MQ)> 11. Our estimates of Zbuiiding for 
passive galaxies {\og{SSFR/Gyr~^)< —2) are shown in Fig. 
1181 the open circle symbolis were derived completely from 
the zCOSMOS dataset, using the GSMF in the first redshift 
bin {ipQ — <j){z ^ 0.24) as reference, and are shown along 
with the region covered using various ETG classification 
criteria. Most of the points have A^>7\/(min, above which 
we are formally complete. Since the local ETG GSMF, de- 
rived by Baldry et al. (2006, 2008) as described previously, 
approaches the observed zCOSMOS GSMF in the first red- 
shift bin for log(A^/A^0)> 10.5 for passive galaxies, below 
this mass the values of Zbuiiding should be considered upper 
limits. We also show Zbuiiding derived by assuming the local 
GSMF for ETGs, along with the region covered using the 
different fractions of ETGs in the 4 central environment 
bins in Baldry et al. (2006). In addition, systematic uncer- 
tainties caused by the different classification methods could 
affect our estimates of Zbuiiding, as for the zCOSMOS data. 

First of all, we note that the "downsizing" trend with 
mass in terms of ETG Zbuiiding is clearly visible. Our data 
are in very good agreement with the values derived by 
Cimatti et al. (2006), who used the de-evolved luminos- 
ity function for galaxies on the red-sequence. The values of 
•^building are in quite close agreement with the end of the 
star formation (80% of the star formed) histories estimated 
by Thomas et al. (2009) for early-type galaxies as a func- 
tion of stellar mass on the basis of their spectral properties 
(T09, in the right panel of Fig. [T8|) . This downsizing trend 
with stellar mass conflicts instead with the prediction of 
semi-analytical models for Zassombiy (e.g., by De Lucia et 
al. 2006, DL06, reproduced here in the right panel of Fig. 
[TS|) . Given the increase with mass of the redshift at which 
the stars are predicted to form in this model (zform, see Fig. 
[18] right panel), the decreasing trend with mass of ^assembly 
is caused by the increasing relevance of merging processes 
(in particular minor merging) with increasing mass (Wang 
& Kauffmann 2008; Stewart et al. 2009; Parry et al. 2009) 
in the SAMs. 

To perform a more accurate comparison with semi- 
analytical models (SAM), we used Zbuiiding derived accord- 
ing to exactly the same definition applied to the data. We 
adopt the GSMF predictions for the three prescriptions 
used in Fontanot et al. (2009): the most recent implemen- 
tation of the Munich SAM (De Lucia & Blaizot 2007, Wang 
et al. 2008a, hereafter WDL08), the MORGANA model 
(Monaco et al. 2007, Lo Faro et al. 2009), and the fidu- 
cial model presented by Somerville (2008, hereafter S08). 
We used the same classification criteria for passive galax- 
ies (log(5S'FiZ/Gyr-i)< -2) as used in the models. The 
predicted GSMFs for each model are shown in Fig. 6 of 
Fontanot et al. (2009). In Fig. [151 we compare Zbuiiding de- 
rived from GSMF SAM predictions (provided by Fontanot) 
and zCOSMOS data using exactly the same definition. For 
each model, we show the region covered by Zbuiiding using 
the original SAM GSMFs and those that take into account 
the 0.25 dex of uncertainties in the mass derivation. In con- 
trast to the data, the semi- analytical models show a reverse 



trend with mass of Zbuiiding, which is similar to the trend of 
•Zassombiy Rcsidual Small differences, for example between 
the determination of mass or SFR for the data and that of 
the models (e.g., the presence of secondary bursts in the 
model SFH) or in the ETG classification criteria are un- 
likely to reconcile the opposite trend with mass. 

Taken at face value, the coincidence of Zform estimated 
by Thomas et al. (2009) and Zbuiiding argues against the 
dominant contribution of mergers (major-|-minor) at z < 1 
in building up ETG galaxies, as instead predicted by semi- 
analytical models (Wang & Kauffmann 2008, Stewart et 
al. 2009). We therefore conclude that the build-up of ETG 
galaxies follows the same "downsizing" trend in mass as the 
formation of their stars, which disagrees with the "upsizing" 
trend predicted in SAMs. 

8.4. Timescales for the quenching of the star formation and 
morphological transformation 

The evolution with redshift in the galaxy number den- 
sity in different mass ranges for three different samples of 
ETG (red {VT =1), spheroidal {ZT =1,2.0), and quiescent 
{log{SSFR/GyT-^)< -1)) galaxies is shown in Fig.M We 
find that they all show a very similar trend (but a different 
normalization), there being an increase in the number den- 
sities of intermediate-mass ETGs with cosmic time from 
z = 1 to the local Universe. This trend is steep for low-to- 
intermediate mass galaxies and is flat and almost negligible 
only for log(A^/A^0)> 11 (< 0.1 dex between z = 0.85 and 
z = 0.25 for log(7W/7Wo)=ll-11.5). The similar evolutions 
of red and spheroidal galaxies was already noted by Arnouts 
et al. (2007), who commented: "If not by chance, this co- 
incidence could suggest that the build-up of the quiescent 
sequence is closely followed or preceded by a morphological 
transformation." In principle, if all the various classification 
methods here adopted are "perfect" , that the number den- 
sity of SED-ETGs is higher than that of the morpho-ETGs, 
at a given mass, suggests that the colour transformation 
(from blue to red) precedes or has a shorter timescale than 
the morphological transformation. From Fig. [THl the delay 
time between the colour and the morphological transfor- 
mation could be estimated, at a fixed pN, to be about 1-2 
Gyr (as shown by the arrow for log(A^/A^Q)=10.7-ll). In 
addition, since the number density of the low- S SFR popu- 
lation is higher than that of the red galaxies (SED-ETGs), 
we estimate that the time elapsed for quiescent galaxies to 
completely switch-off their SFR and become red is of the 
same order (about 1-2 Gyr, shown as an arrow in Fig. I19p . 
Wolf et al. (2008) reach a similar conclusion in the cluster 
A901/2 at z ^ 0.17: the rich red-spiral population at inter- 
mediate mass {\og{A4/ A4q)= [10,11]) is more accurately 
explained if quenching is a slow process and morphological 
transformation is delayed even more (see also Skibba et al. 
2009). The small contribution to the assembly of the red- 
sequence by the class of post-starburst galaxies (most of 
which are assumed to quench their star formation rapidly) 
provides additional support to a scenario of delayed/slow 
quenching (Vergani et al. 2009). 

However, we are aware that the colour-selected red 
galaxies may be contaminated by dusty starbursts. We re- 
call here that, given our 24 /im flux limit of ~ 0.3 mJy, only 
5% of red galaxies (SED-ETGs) were detected at 24 ^m 
(all having log{SSFR/GyT-^)> -2), and about 23% have 
—2 <log(S'5'FiZ/Gyr^^)<-l, i.e., they are not completely 
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Fig. 19. Number density evolution for ETGs in different 
mass ranges: red galaxies (red filled triangles for VT—1), 
spheroidals (dark magenta empty circles for ZT— 1+2.0) 
and low-SSFR (magenta filled circles for \og{SSFR)< -f). 
Black filled squares refer to tfie global population. Lines 
have the same meaning as in Fig. 1111 In the second panel, 
we also indicate by 2 arrows the evolutionary time sequence 
from quiescent to red to spheroidal galaxies (see text). 



"dead". Bundy et al. 2009, through MIPS stacking 
analysis, show that the implied SFRs of red galaxies 
are on average at least an order of magnitude be- 
low the star-forming population. On the other hand, 
there is a non-ncgiigible blue population (^ 20%) among 
elliptical galaxies, for which therefore the morphological 
transformation could precede or is faster than the colour 
transformation (see below for a discussion). Therefore, we 
also checked that the number density of red passive galax- 
ies (VT =1 and \og{SSFR/GyT-'^)< -2), i.e., those that are 
"red & dead" , is higher at any given mass and redshift than 
the number density of passive elliptical galaxies (ZT=1,2.0 
and log(S'5'Fi?/Gyr~-^)< —2), and that at a given number 
density and mass they are delayed by 1-2 Gyr. We conclude 
that, at least statistically, the morphological transformation 
takes a longer time than for a given galaxy to become "red 
and dead" . 

But which physical processes induce these transforma- 
tions? All the properties SSFR, colours, and photomet- 
ric types depend on the star-formation history, while mor- 
phological types mainly reflect the dynamical history of a 
galaxy. The transformation from active blue galaxies to red 
passive galaxies, as well as the transformation from spi- 
ral disc-dominated to bulge-dominated systems, could be 
driven by a variety of processes such as the ageing of their 
stellar population or the fading of the disc, because of the 
exhaustion of its gas reservoir, gas stripping, or the quench- 
ing of the SFHs, for example by AGN feedback (Menci et al. 
2006; Monaco et al. 2007; Bower et al. 2006), or by truncat- 
ing the gas accretion from infall or from satellite galaxies. 
In addition, "starvation" (or "strangulation") is expected 
to affect a galaxy star formation history on a quite long 
timescale and therefore to cause a slowly declining activity: 
when a galaxy is accreted onto a larger structure, the gas 
supply can no longer be replenished by cooling, which is 
suppressed (Larson et al. 1980). This process is an impor- 



tant element of semi-analytic models of galaxy formation. 
Some of the aforementioned processes are also obviously 
environment-dependent . 

On the other hand, the strong increase in the num- 
ber density of intermediate-mass pure morphological ellipti- 
cals (ZT=1) or "bona- fide ETGs" (red passive spheroidals) 
provides support for dynamical processes, such as the merg- 
ing between two low-mass galaxies (Toomre & Toomre 
1972, Barnes 1992). Mergers are included in standard semi- 
analytic models of galaxy formation and represent indeed 
the main channel of formation for bulges. Menci et al. 
(2002) also included in their hierarchical semi-analytical 
models binary aggregations of satellite galaxies inside a 
common halo. However, if merger events were responsible 
for morphological transformations, we would expect to ob- 
serve a rapid evolution from blue to red after the forma- 
tion of the bulge/elliptical, i.e., after the burst triggered 
by the merger (Mihos & Hernquist 1996). We find, instead, 
that the colour transformation precedes the morphological 
one, and that only a population (~ 20%) of massive blue 
spheroidal galaxies exists, which could be the end-product 
of merger events. The massive tail of this population has 
some signature of an AGN component in its stacked spec- 
trum (see Fig. [77]) , which may be responsible for quenching 
the SFR, allowing their transition to red ellipticals with 
cosmic time. The study of this particular population is be- 
yond the aim of this paper and will be the subject of future 
zCOSMOS papers. 

Our results, therefore, suggest that dynamical pro- 
cesses, other than mergers, are the main mechanisms re- 
sponsible for the morphological transformation that fol- 
lows the colour transformation. For example, bulge growth 
through disc instability could be efficient (CaroUo et al. 
2001, Bower et al. 2006, Dekel et al. 2009). Disc instabil- 
ity is also a possible formation channel for bulges in SAMs 
(De Lucia et al. 2006, Bower et al. 2006, Parry et al. 2009). 
Dynamical observations of massive galaxies at z ~ 2 with 
SINFONI (Genzel et al. 2006, 2008) also suggest that disc 
instability can produce massive bulges, even in the absence 
of major mergers. Indeed, dynamical friction and viscous 
processes at z ~ 2 proceed on a timescale of < 1 Gyr, which 
is at least one order of magnitude faster than in z ^ disc 
galaxies. In overdense environments, such as clusters, ram 
pressure may also play a role, demolishing the star-forming 
disc (Cortese et al. 2006). Recently, Bundy et al. (2009) 
demonstrated that the colour and morphological transfor- 
mations should proceed through several separate stages and 
explored the strengths and weaknesses of several more so- 
phisticated explanations, including environmental effects, 
internal stabilization, and disc regrowth by means of gas- 
rich mergers. 

9. Comparison with models 

To establish the cause of the "age-downsizing" of ellipti- 
cal galaxies, according to which the stars in more mas- 
sive galaxies formed earlier and over a shorter period than 
those in less massive galaxies, Cattaneo et al. (2008) dis- 
cuss a model in which star formation shutdowns in dark 
matter haloes of mass above a critical value of ^ 10^^ -Mq- 
However, this model predicts an evolution in the massive 
tail of the GSMF of red galaxies that is more rapid than 
that of intermediate-mass red galaxies from z — 1 to z — 0, 
which is in contrast, at least qualitatively, with our results 
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Fig. 20. Comparison between observed GSMFs (lines and points) with the regions covered by the predictions (shaded 
colour regions) of three semi- analytical models (Morgana, WDL08, and SOS from Fontanot et al. 2009), at increasing 
redshifts (from top to bottom panels). Left panels: original GSMFs from SAM. Right panels: SAM GSMFs take into 
account 0.25 dex of uncertainties in the Ai derivation. 



(see Ilbert et al. 2010 for a comparison between model and 
data). 

The same behaviour is clearly visible in Fig. 6 of 
Fontanot et al. (2009), which shows that the evolution of 
the GSMF for passive galaxies with redshift follows ex- 
actly the opposite trends with mass and cosmic time ob- 
served in the data (compare with our Fig. I16p . They also 
conclude that a robust prediction of SAM models seems 
to be that the evolution of less massive galaxies is slower 
than for more massive ones, i.e., the models do not pre- 
dict "mass-downsizing" but rather the opposite behaviour 
(sometimes called "upsizing" ) . Here we recall that the three 
semi-analytical model renditions considered in Fontanot 
et al. (2009) are indeed unable to reproduce the "mass- 
assembly downsizing" evolution of observed ETG GSMF, 
as parametrized for example by the median building red- 
shift described in the previous section (Fig. [T8|) . 

In Fig. [201 we also show the direct comparison of zCOS- 
MOS ETG GSMF with the SAM model predictions for 
passive galaxies (log(S'S'Fi?/Cyr~-^)< —1). Hatched region 
show the range covered by the three models from Fontanot 
et al. (2009) for the original GSMFs (left panel) and con- 
volved with 0.25 dex uncertainties in the mass estimate 
(right panel). A more detailed and quantitative compar- 
ison between each model and data is beyond the aim of 
this paper. We first note that the shape of the predicted 
GSMFs differs from that of the observed GSMFs in all the 



mass and redshift ranges explored. In particular, the mod- 
els are unable to recover the steep decline with decreasing 
mass in the number densities of passive ETGs, but instead 
overproduce their number densities increasingly with de- 
creasing mass for log{M/MQ)< 10.3. This is because in 
the models, intermediate-mass galaxies form too early and 
are too passive and red at late times. In contrast, at the 
massive mass end (log(A4/A^o)> 11) the original mod- 
els have the well-known problem of underestimating the 
number density of massive galaxies at high redshift (for 
z > 0.6 in zCOSMOS dataset). We show that the convo- 
lution with stellar mass uncertainties (0.25 dex) only par- 
tially solve the problem at high redshift and high masses 
{\og{M/M(7))> 11.3), while models continue to underpre- 
dict GSMFs around log(7W/A^0)- 11 - 11.3, where they 
fail to reproduce the GSMF shape, and overproduce mas- 
sive galaxies at low redshifts, as already pointed out by 
Fontanot et al. (2009). To qualitatively reconcile data and 
models at the massive-end, the mass uncertainties would 
have to increase rapidly with redshift, which is not the case 
for most current surveys with accurate and multi-band pho- 
tometry that extends to near-IR rest-frame wavelengths. 

We conclude that the three semi-analytical model ex- 
plored here (WDL08, MORGANA, and SOS) are unable 
to completely reproduce either the observed shape of the 
GSMFs of passive galaxies or the evolutionary trend. In 
particular, SAMs overpredict the low-mass ETG number 
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densities and, even if the uncertainties in A4 determination 
partially assount for the underprediction at high masses and 
high redshifts, they predict an evolution with cosmic time 
that is of an opposite trend with mass to that observed for 
the data, i.e., they arc unable to reproduce the later forma- 
tion and more rapid evolution in number densities of the 
lower mass galaxies ("mass-assembly downsizing"). 

10. Summary and conclusions 

We have investigated the evolution of the galaxy stellar 
mass function up to ^; = 1 in the zCOSMOS 10k bright 
spectroscopic survey using VIMOS spectroscopy 8500 
galaxies with 15.0 < I < 22.5 over 1.4 deg^) and multiband 
photometry (from UV to near-IR). Our main observational 
results can be summarized as follows: 

— Unveiling bimodality in the GSMF: The shape of the 
GSMF is more accurately reproduced by two Schechter 
functions at least up to 2; ~ 0.55, the maximum redshift 
at which we can explore the GSMF upturn due to our 
mass limit. This shape is linked to the bimodality in 
the galaxy properties of the ETG and LTG populations, 
which dominate the GSMF, respectively, at high masses 
{\og{M/MQ)> 10.5) and at low/intermediate-masses 
up to the highest redshift explored (z ~ 1). 

— Mass-assembly downsizing: We find a continuous, rapid 
increase with cosmic time from z ~ 1 to 2; = in 
the global GSMF for log{M/MQ)< 11 (by a factor 
of ~ 2 for log(A^/A^0)^ 10.5) and a more slow in- 
crease (< 15%) for log(A^/Al0)> 11. A similar but 
even stronger trend is inferred for the ETG population, 
i.e., massive galaxies build-up their stars and mass ear- 
lier than lower mass galaxies. We show that the Zbuiiding 
of ETGs increases rapidly with Ai from Zbuiiding^ 0.2 
to ~ 0.8 for log(A^/A^0)= 10, 10.8, respectively, and is 
significantly higher {z ^ 1) for massive ETGs. 

— "SFH-evolved GSMF" and mergers: Assuming only a 
growth in stellar mass driven by exponentially decreas- 
ing star-formation rates, we find that the predicted 
"SFH-evolved GSMFs" agree with the observed GSMFs 
with differences of at most 20-40%, even if they predict 
redder colours than observed, in particular at low red- 
shift. Therefore, mass assembly is controlled primarily 
by SFH, which can explain most of the evolution of the 
GSMF at intermediate-low M {log{M/MQ)< 10.6). 
The low residual evolution of the GSMF is consistent 
with 0.16 merger per galaxy per Gyr, on average, with 
a hint of a decrease with cosmic time but no clear de- 
pendence on mass. Major merging events contribute in 
terms of fewer than 0.1 merger /gal/ Gyr (~ 0.06 on av- 
erage) and are, if not marginal, not the dominant evo- 
lutionary process. 

— A^cross and the evolution of the LTGs: The intersection 
(Atcross) between the GSMFs of the ETG and LTG 
populations decreases with cosmic time. This is due 
mainly to a clear increase with cosmic time in the num- 
ber density of ETGs, regardless of their classification 
method, for log{M/MQ)< 11. In contrast, for LTGs 
the number density of blue or disc-|-irregular galaxies 
shows only a mild or negligible evolution, while the most 
extreme population of star-forming galaxies (those with 
high specific star formation) decreases rapidly in num- 



ber density with cosmic time, in particular at high mass 
(log(7W/7W0)> 10.3). 

— Flow to build-up the red population: We have quanti- 
fied the fiow rate from bhie to red population at dif- 
ferent redshifts and masses. We have found a growth 
rate in number and mass density of the red galaxies 
of about a few 10^'' gal/Mpc'^/Gyr/logAI and a few 
10*^ - 10'' 7W0/Mpc3/Gyr/logA4 for \og(M/MQ)< 11, 
and lower in value above this mass. The correspond- 
ing fraction of blue galaxies that, at any given time, 
are transforming into red galaxies per Gyr is on av- 
erage ~ 25% for log{M/MQ)< 11. This fraction and 
the growth rate for red galaxies increases with cosmic 
time at \og{M/MQ)= 10.4 - 10.7 at least for z < 0.7, 
while it seems to decrease with cosmic time at higher 
M {logiM/MQ)= 10.7-11). 

— ETG evolutionary sequence: Using different methods to 
classify ETGs, we have found that the number den- 
sity of quiescent galaxies is higher than that of red 
galaxies, which is even higher than that of spheroidals 
for log{^A/^AQ)< 11. This behaviour suggests that a 
transformation occurs from quiescent to red objects, 
as inferred mainly from SFHs, before or on a shorter 
timescales than the morphological tranformation, which 
is most likely primarily driven by dynamical processes. 

— Comparison with SAM: Semi-analytical models are un- 
able to fully account for GSMF evolution, in particular 
for ETGs. They predict a different shape at any mass 
and redshift explored, there being a strong excess of low- 
mass passive galaxies and a deficit of high-mass galax- 
ies at high-z, disagreement that can be only partially 
resolved by the bias introduced by uncertainties in the 
mass determination. Furthermore, they predict an "up- 
sizing" evolution in number density with cosmic time, 
i.e., mass assemply occurring slower and earlier in low- 
mass galaxies, which is the opposite of the "downsizing" 
trend observed in zCOSMOS data for the GSMF and 
quantified by ^building- 

We attempt to summarise our results within a galaxy 
evolutionary scenario that emerges from several observa- 
tional studies. It is almost established that there is an 
age-downsizing of galaxy stellar populations, according to 
which less massive galaxies contain younger stars, which 
is also applicable to similar spectral/morphological types 
both among ETGs (Nelan et al. 2005; Thomas et al. 2005; 
Graves et al. 2007; Fontana et al. 2004) and LTGs (Noeske 
et al. 2007a,b). Furthermore, we find that downsizing is also 
applicable to the mass-assembly, i.e., massive galaxies as- 
semble their mass in a single object earlier than low-mass 
galaxies, both for the global population and for ETGs, re- 
gardless of their definition (in terms of their colours, mor- 
phology, or SFR activity). 

We therefore propose a scenario where intermediate- 
mass galaxies {log{A4/ A4q)= 10 — 11) decrease their star 
formation activity to intermediate values gradually (with 
an e-folding time r = 1 — 3 Gyr) and finally to that of qui- 
escent galaxies, by means of the exhaustion of either the gas 
reservoir or cold gas accretion, or quenching due to AGN 
feedback. Noeske et al. (2007a, b) found that the e-folding 
time is mass-dependent, being longer for less massive galax- 
ies, which also have a later onset of the SFH. Therefore, if 
no other mechanism intervenes to again switch on the star- 
formation activity (e.g., re-juvenation hypothesis, Hasinger 



L.Pozzetti, et al.: zCOSMOS: Galaxy bimodality in the stellar mass function 



23 



et al. 2008) , they evolve to have red colours on a timescale 
of 1-2 Gyr and finally undergo a dynamical morphological 
transformation into spheroidal galaxies. 

In this scenario, massive galaxies (log(A^/7Wo)> 11), 
in particular ETGs, are already in place at z = 1, and their 
GSMF exhibit negligible evolution thereafter. In contrast, 
the number density of intermediate-mass ETGs continues 
to increase with cosmic time (by a factor ~ 4 — 2 from 
z = 0.7 to z = 0.2 for \og{M/MQ)= 10.3 - 10.7) with 
a growth rate of a few 10"'' gal/Mpc^/Gyr/logA^ and a 
few 10^ A^o/Mpc^/Gyr/logA^. We have found that the 
most likely evolutionary process to explain our results is 
the transformation from blue to red galaxies by means of 
the smooth decrease or the quenching of the SFR (as de- 
scribed above) , followed by a morphological transformation 
probably by means of internal dynamical processes, such 
as disc instabilities. Major merger events are unlikely to 
be able to explain why colour transformation precedes the 
morphological one. In contrast, blue galaxies at low mass 
(log(A^/A^0)< 10) continue to grow in stellar mass, be- 
cause of their high SSFR, and replace intermediate-mass 
LTGs, whose density remains almost constant in time. 
We have found that only a small fraction of low-mass 
(log(A4/A^0)'~ 10) blue galaxies need to transform into 
red galaxies to explain the observed ETG growth rate at 
least at z > 0.4. Merger events (minor and major) con- 
tribute with a rate of ^ 0.16 merger/gal/Gyr on average 
(of which fewer than 0.1 are major) to form galaxies at 
redshifts lower than z ~ 1. 

Following a schematic diagram, as originally proposed 
by Faber et al. (2007), this scenario is presented in Fig. 
\n\ The dominant processes are the transformation with 
cosmic time of blue star-forming galaxies of progressively 
lower mass to red passive galaxies according to their SFH. 
Different evolutionary tracks are plotted between 0.1 and 
10 Gyr and refer to constraints on their inferred SFHs 
with different timescales (longer for less massive galaxies) 
and normalized to the total final mass (Mtot)- As exam- 
ples that only qualitatively reproduce our data, we show 
tracks for models with (r; log(7Wtot/A^o))= (0.1; 11.5), 
(0.3; 11.0) [solid red curves]; (0.6; 10.7), (1;10.3), (2; 10.0) 
[short dashed green curves]; (5; 11), (10; 10), (15; 9) [long 
dashed blue curves]. 

We conclude that the build-up of galaxies, and ETGs 
in particular, follows the same "downsizing" trend with 
mass (i.e., occurring earlier in high-mass galaxies) of the 
formation of their stars, which is the opposite of that pre- 
dicted by current SAM models. The SAM LCDM mod- 
els, with rapid galaxy growth, possibly overestimate 
the merger rate within hales (and underestimate 
the time-scale for merging). Our analysis also sug- 
gests that morphological transformation should oc- 
cur on longer timescales than colour transforma- 
tion. Merging processes are not marginal but their fre- 
quency (~ 0.16 merger/gal/Gyr of which < 0.1 are major) 
appears to decrease with redshift, while no strong depen- 
dence on mass has been found, in disagreement with SAM 
models (Wang & Kauffmann 2008, Stewart et al. 2009). 
The prevailing galaxy formation and evolutionary scenario 
is therefore moving towards one in which a smoother evolu- 
tion in mass growth and star formation (due to accretion of 
cold gas) and minor merger play a more important role than 
major merging events. In this scenario, however, it appears 
quite difficult to explain the "fine-tuning" of the transfor- 



red sequence 




Fig. 21. Schematic scenario for galaxy evolution. Different 
ellipses refer to different redshifts and galaxy populations, 
while different evolutionary tracks are plotted between 0.1 
and 10 Gyr and refer to SFHs with different timescales and 
normalized to the total final mass ((r; log(A^tot/A^o))= 
(0.1; 11.5), (0.3; 11.0) [solid red curves]; (0.6; 10.7), (1; 10.3), 
(2; 10.0) [short dashed green curves]; (5; 11), (10; 10), (15; 9) 
[long dashed blue curves]. 

mation of low-mass into intermediate-mass blue galaxies re- 
quired to ensure that their number density remains almost 
constant at any given mass. In this scenario, we expect 
negligible evolution of the galaxy baryonic mass function 
[GBMF) for the global population at all masses and a de- 
crease with cosmic time in the GBMF for the blue galaxy 
population at intermediate-to-high masses. 
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